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I. INTRODUCTION 



Every object in our environment emits thermal radiation. The peak energy of this 
radiation is characteristic of the object's temperature and. for temperatures coimnonly 
found in our environment, is concentrated mostly in the infrared region of the 
electromagnetic spectrum. This physical phenomenon has found many applications 
since the second world war as it enables passive detection of targets by their self- 
emission alone. 

Passive search, detection, identification and tracking have always been of great 
importance in naval warfare. Also, the necessity to improve our combative capabilities 
in areas where radar performance is degraded (low Hying small cross section targets) is 
imminent. Infrared technology promises much in this area. The rapid growth of 
technology in solid state integrated circuits, microprocessors, infrared detectors and the 
better understanding of infrared physics during the recent decades has made it possible 
for sophisticated infrared systems to be designed and developed. 

A very important step in designing an infrared optical system is the knowledge of 
the spatial and spectral radiances of the target and its background as well as the 
atmospheric elTects on the emitted thermal radiation. Almost always, targets appear 
against a background that complicates the detection process since radiation from 
objects surrounding the target is confused with radiation emitted or reflected from the 
target. In the analysis of thermal imaging systems the targct-to-background radiance 
contrast is frequently described by the equivalent temperature dilTerence. In the case of 
observation at low elevation angles above the sea surhice the surface radiance is 
strongly influenced by the reflection of the radiation emitted by the sky. The actual 
thermometric temperature difference is no longer a good representation of the radiance 
contrast. An "effective temperature difference" can however be defined taking into 
account the reflected radiation. Any new information and additional data about 
background scenes will improve the understanding of these elTects and may be used to 
optimize computer algorithms or to implement new or improved techniques in signal 
processing to minimize background effects, thus improving system capabilities and 
optimizing performance. 
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This thesis describes the collection and analysis of thermal images of an 
instrumented ship target and its background under various well defined environmental 
conditions, followed by the development of a model for the effective target-to- 
background temperature difference. The target was the oceanographic ship R V 
"POINT SUR" and the background was sea surface and sky. The data were collected, 
stored, and processed using the "AGA 7S0 THERMOVISION", a liquid nitrogen 
cooled thermographic device with digital recording and data analysis hardware. .All 
measurements were taken from a distance of about 6o0 meters using the Mercury 
Cadmium Telluride (MgCdTe) detector of the Thermovision 7S0 scanning system in the 
8-14 pm spectral region. Five experiments were conducted during an operational 
oceanography student cruise on mid May 1986 and an adequate number of thermal 
image pictures was taken for various environmental conditions on certain dates. 

.Analysing these thermal images using the measured ship surface emissivity and 
the AGA software package, both the radiance and the effective temperature 
distributions of the ship were obtained. The atmospheric transmittance for this purpose 
was computed using the LOWTR_AN 6 code, and the emissivity was measured from 
samples of the ship paint. The effective sea surface temperature was also obtained from 
these pictures. It was observed that the apparent sea surface temperature was affected 
by reflected sky radiation, and a model was developed to describe the effective target- 
to-background temperature difference under varying environmental conditions. 

Later, three more experiments were conducted in mid November 1986 in order to 
obtain further, more accurate, data under different environmental conditions, after the 
Thermovision had been calibrated. This was done to check the previous results and 
improve the information on the reflected sky thermal radiation and the effective target- 
to-background temperature difference relation. 

The present work consists of eight chapters. In the first chapter we state the 
background problem. In the second chapter we discuss the background radiances which 
affected our measurements, the marine background and the sky background. The 
atmospheric propagation problem and the LOWTRAN code used to find the 
atmospheric transmittance for the various environmental conditions are discussed in 
chapter three. Chapter four provides a familiarization with the equipment used to 
collect the data. In the fifth chapter we describe the calibration problem and discuss 
the derivation of the empirical calibration relations. In the sixth chapter we discuss 
how the thermal data we had collected were analyzed to give the temperature and 
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radiance distribution of the target. A further analysis of these data, to see how the sea 
background radiance was alTected by reflected sky thermal radiation, is given in 
Chapter seven. In this chapter we also calculate the apparent sea surface temperature 
based on environmental parameters and compare the target-to-background temperature 
difference sensed by our detector with the actual thermometric temperature difference 
measurements. In the last chapter we discuss our conclusions and recoinendations for 
further work. The discussion and analysis of the data obtained during the mid 
November measurements are given in Appendix A. In Appendix B the TI-59 program 
used to get the in-band radiant flux for the radiance calculations is explained. 

The results of this work will lead to a better understanding of the background 
effects and may add some information to be used for the implementation of techniques 
in signal processing. 
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II. TARGET AND BACKGROUND 



A. BACKGROUND RADIATION NOISE 

Background radiation noise is noise generated in an optical detector by photon 
flux reaching it from all points other than the target. This kind of noise usually is 
separated into two components [Ref 1]: 

1. Temporal, which is due to fluctuations in the background einittcd photons. 

2. Spatial, which is due to background structural details. 

The temporal background noise is generated by the fluctuation in the number of 
incident photons at the detector. This noise usually is referred to as background 
photon noise and is not directly related to the target. Detectors sensitive enough to 
experience background photon fluctuations noise are called BLIP ( Background 
Limited Photon Detectors ). 

In this work we are interested in the spatial component of the background noise, 
produced by the structural content of the background and mainly related to the target. 
Therefore, from now on whenever we refer to background noise we will mean the 
spatial component. 

Since a target is always surrounded by some sort of background it is obvious that 
radiation reaching our detector from this background is confused with the radiation 
emitted or reflected from the target. This is a serious problem that complicates the 
detection process and must be well understood for designing a high pcrforance optical 
or infrared system. [Ref 1] 

Background radiation is due to thermal emission and reflected or scattered 
incident radiation. Many parameters afTect this kind of radiation. Of high importance 
among them are the changes in emissivity and reflectance of the various background 
surfaces under varied environmental conditions. All the above make the background 
problem complicated and the modeling of a universal background function impossible. 

Background radiance can be considered mainly as two categories: 

1. Terrestrial background 

2. Celestial background 

Terrestrial background refers to radiance seen by infrared systems working within 
the atmosphere of the earth. This kind of background consists of radiation that several 
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surfaces of the earth are emitting or reflecting { oceans, clouds, vegetation, etc. ). It 
also includes sunlight scattered by particles in the atmosphere and atmospheric- 
molecular emission ( sky radiance ). 

Celestial refers to that background seen by a system working in space and 
consists of radiation from the stars and the sun reflectance from planets and space 
debris. 

Of particular interest in this work are the background radiances of the sea surface 
and the sky. Detailed discussion of these backgrounds and their spectral radiance 
characteristics, especially for the spectral region of 8 to 14 ^m. where we have taken 
our measurements, is given in Chapter III. 

Since target and background radiances travel through the atmosphere to reach 
the detector, they are heavily influenced by the atmospheric absorption under various 
meteorological conditions. Detailed analysis of the atmospheric effects on thermal 
radiation is given in Chapter IV. 



B. BACKGROUND RADIATION CONTRAST 

The target-to-background temperature difference is an important factor in the 
target detection probability of an infrared system. The higher the temperature 
difference between a target and ' its background, the greater the distance and the 
probability of detection. 

This temperature difference or radiance difference is usually e.xpressed with the 
quantity Radiation Contrast. Cj^. defined in [Ref 2] as 



■ 



\Vt - Wg 



\Vj + Wg 



where W-g. and Wg are the target (T) and background (B) radiant emittances (W m“). 

The 8 to 14 pm wavelength radiation contrast versus target-to-background 
temperature difference for different background temperatures is shown in Fig 2. 1 . 

When the temperature difference between the target and its background is very 
small, or goes to zero, the radiation contrast goes to zero and the detection becomes 
ver>- difficult if not impossible. This effect is known as the washout effect and happens 
for "ambient temperature" targets twice a day due to the solar cycle. 
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Figure 2.1 Radiation contrast for tlie S-14 um region. 



In [Ref. 3] Hudson describes an experiment where the importance of the contrast 
effect is shown. In that experiment a vehicle parked in an open field was viewed by an 
infrared optical system over a period of 24 hours. During that observation, contrast 
variations were found. In the afternoon the contrast was positive, since the vehicle 
heated by the sun was warmer than its background. Later in the early evening hours 
the contrast was higher, since the vehicle with larger thermal capacity cooled more 
slowly than the background. .A.S the night progressed, the contrast passed through zero 
and went negative, since the vehicle continued to cool rapidly and at some time became 
colder than its background. During the morning the negative contrast was still present 
because the background was warnting more rapidly. .After some hours the heating of 
the vehicle was enough to give a period of zero contrast and then turned to positive 
values again. 

Hence, twice in a 24 hour period we are lacking sufficient contrast level between 
the target and its background, which means that the detection of any target at these 
periods is not possible, and since this effect can not be eliminated we have to utilize 
other methods of detection. 

C. BACKGROUND SUPPRESSION 

Target and background have different spectral radiances and different geometries, 
differences that make it easier for them to be distinguished. Whenever a spectral 
radiance difference exists between a target and its background an optical spectral filter 
can be used to reduce the unwanted background signals. Even if spectral, filtering 
increases the target to background ratio, this is rarely sufficient to render a system 
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operational since usually the irradiance of the background is much higher than that of 
the target. A common characteristic of targets is that they are always smaller in 
angular extent than their backgrounds, a typical example being a ship against the sea. 
The difierences in geometry between targets and background can be expressed as 
differences between their spatial frequency components. For such geometry dilTerences. 
spatial filtering provides a good suppression of backgrounds of this extent relative to 
the target radiance. 

Spatial filters include all types of reticles discussed in detail in [Ref 3]. Spatial 
filtering is supplemented by spectral filters that define an effective spectral bandpass 
for the system and electrical filters that work on an average background threshold level 
to chop the background dc signal. 

It. is easily understood that spatial filtering is the basic suppression scheme 
mainly giving discrimination of the target in an extended background. But for many 
military’ systems, depending on the tasks they perform ( identification, tracking, etc.), 
this is not enough. We have interference from other smaller kinds of background such 
as unwanted objects in the search area, for example buildings, rocks, clouds, etc. or 
unexpected objects crossing our detector's field of view such as birds, different kinds of 
ships, etc. 

Therefore, more sophisticated and complicated spatial filtering systems have to 
be used for all the different purposes. These systems are supported by improved 
electronics and computer software, using different implementing techniques to minimize 
the background effects and optimize the system's performance. 
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III. BACKGROUND RADIANCE 



A. MARINE BACKGROUND 
1. Introduction 

The sea surface infrared radiance (Lj) is the sum of the radiance caused by its 
thermal emission (L^^.(l-R)) from below the surface and the radiance due to reflected 
incident radiation (Lj.R) (R: reflectance of the sea surface). 




Lj = L,,.(l-R) + Lj.R 
Figure 3.1 Sea surface radiance. 

Major factors that determine the marine background character are the 
following [Ref 4]: 

a. Infrared optical properties of sea water 

b. Geometry of the sea-surface and wave slope distribution 

c. Sea water temperature distribution 

d. Properties of the bottom material . 

2. Infrared Optical Properties of Sea Water 

Thermal radiation from the sea surface occurs at wavelengths having their 
maxima close to the 10 pm region. .‘Xs we can see in Fig 3.2 taken from [Ref oj. at 
these wavelengths water is essentially opaejue in thickness of about U.02 mm. .At 
wavelengths above 3.5 pm only the upper layers of the sea, of thickness about 0.01 cm 
and less, determine the sea surface thermal emission. 
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Figure 3.2 Sea water absorption for selected wavelength regions. 

The radiation emitted from the sea surface is dependent on the underwater 
temperature profile, and is little affected by surface contamination. Since the emissivity 
is equal to the absorptivity the sea surface radiance is dependent on the distribution of 
the absorption coefficient K of the sea water as shown in Fig 3.3 . 



Wavelength 




Figure 3.3 Absorption coeflicient K of sea water versus wavelength. 
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Subsurface scattering of incident radiation is not significant in the infrared 
range, though it happens in the solar range (0.35 - 3.0 }im). In that region the 
absorption coefTicient has a ver>' low value and radiation penetrates deeply into the 
water and is scattered. 

Reflected incident radiation is composed of solar radiation and sky thermal 
radiation . Solar radiation dominates at short wavelength and of course is present only 
in the daytime. In the region from 2 nm to -! pm reflected solar emission especially at 
small angles of incidence is greater than the thermal emission of the sea surface. In 
comparison to this region only about O.OS percent of the sun radiant energy is found in 
the region of S [im to 1-1 }im wavelengths, which means that only reflected sky thermal 
radiation influences this region by day and night. 

The influence of reflected sky radiation was recognized by Me Alistcr [Ref 6] 
who found, for a clear sky at night, an apparent sea surface temperature a few tenths 
of a °C colder than the actual temperature. 

Reflectance, transmittance, emissivity and index of refraction for water, arc 
shown in Fig 3.4 to Fig 3.7 [Ref 4]. Mere wc have to mention that distilled water and 
sea water transmittance do not have any significant difference. 

For the purpose of this thesis, where wc have worked in the S to 14 jim 
wavelength region, we have taken into account only reflected sky thermal radiation. 




Fisure 3.4 Transmittance and reflectivitv of 0.003 cm of sea water. 



HEKLECTrVlTY AND EMISSIVITY (U 




Figure 3.5 Index of refraction of water. 




Figure 3.6 Reflection from water surface at different incidence angles. 




Figure 3.7 Reflectance and emissivity of water. 
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3. Geometry of the Sea Surface 

The slope distribution of the waves is an important factor in the reflection and 
refraction of incident radiation on a roughened sea surface at various wind speeds. In 
Fig 3.S, taken from [Ref. 7], is shown the reflectance of solar radiation from a flat sea 
surface ( ff = 0 ) and from seasurface roughened by a Beaufort 4 wind (a = 0.2 ), 
where c is the square root of the variance of the wave slope distribution. In this figure 
we can see that for an average rough sea ( wind speed of Beaufort 4 ) the reflectance of 
the sea near the horizon goes to about 20^o where as the reflectance of the flat surface 
is about 90% and above. Consequently we can see how high the emissivity ( 1 - R ) of 
the roughened sea surface is near the horizon. Only at higher angles of incidence, about 
40”, does the reflectance of the roughened sea surface go slightly higher than the 
reflectance of the flat surface. 




Figure 3.S Reflectance from a flat ( a = 0 ) and roughened sea surface (cr = 0.2 ). 
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In [Ref. 7J Charles Cox and Walter Munk studied this phenomenon for solar ■ 
radiation in the visible range. They assumed unpolarized sun light reflected once no 
double scattering. They used the index of refraction for uncontaminated sea water as n 
= 1.33S. for calculating the reflection coefTicient. Information is lacking for the 
infrared wavelength region and for night conditions, but since the sky radiance 
variation is almost the same during the day and night, and the reflectance of the sea 
water in the visible region of electromagnetic spectrum at zero angle of incidence is 
almost equal to the average reflectance in the infrared region from 2 pm to 15 pm ( see 
Fig 3.4 ) these curves are considered helpful. For the purpose of this thesis, we have 
used these curves shown in Fig 3.8 from [Ref 7] after making a correction for the index 
of refraction in the infrared range ( 8 to 14 pm ) we have used in our work. .A, detailed 
discussion of this is given in Chapter VIII. 

The problem with these curves is that for the dashed portions near the 
horizon, above 80® angle of incidence, the function being used is not strictly applicable 
because of shadowing and multiple reflections. Also the lower and upper branches of 
the ( (7 = 0.2 ) curve represent two assumptions used in [Ref 7] for the multiple 
reflection efTect. The basic function used works for light reflected only once as we 
mentioned before. True values are expected according to C.Cox and W..Munk to lie 
between the indicated limits. 

The curves in Fig 3.9 show the sea surface radiance N(p) divided by the sky 
radiance Nj(0) at the zenith, due to sun light scattering by air molecules, as a function 
of the vertical angle p. for a clear sky and an overcast sky, and is also taken from 
[Ref 7]. 

4. Sea Water Temperature Distribution 

The sea surface temperature determines the emitted thermal radiation of the 
sea. This temperature is about 0°C (mainly about d^C) in arctic regions and increases 
reaching about 29°C near the equator. Currents may produce anomalies of some 
degree in temperature. However the thermal emission of the sea surface has its 
maximum peak near the 10 pm wavelength. .As we have said before in the infrared 
optical properties section, the temperature of the upper layers of the sea surface, of 
thickness about 0.01 cm. characterizes the thermal emittance. 

Under evaporation conditions the temperature of the upper layer of 0.1 rmii of 
the sea surface was measured [Ref 4] as colder than the water temperature a few 
centimeters below the surface of about 0.6°C. In Fig 3.10 [Ref 4] are shown 
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Figure 3.9 Sea surface radiance divided by sky radiance at zenith. 

measurements of typical conditions showing the sharpest gradient in the upper layer of 
1 mm. 

5. Properties of the Bottom Materials 

These properties are important in regions of short wavelengths where 
radiation can penetrate deep into the water causing subsurface scattering. Such 
properties are not important for the work done in this thesis, since we have worked in 
the region of longer wavelength ( S ].im to 14 fim ). 



26 



TEMPERATUaE (°C) 
12 13 




12 13 

TEMPERATURE (®C) 



Figure 3.10 Thermal structure of sea boundar>' layer. 

B. SPECTR.\L RAD1.A..NCE OF THE SKY 
1. Introduction 

Sky background radiance in the infrared is caused by two dilFerent factors, the 
scattering of radiation from the sun and the thermal emission of the atmospheric 
constituents. The scattering factor is important in the short wavelength range, close to 
the visible, and is present only in daytime. The thermal emission happens at 
wavelengths longer than 4 jtm and is present day and night. 

.An idealized form of this contribution to the sky spectral radiance is given in 
Fig 3.11 [Ref. S]. This figure shows the spectral radiance of a blackbody at 300° K, 
close to the temperature of the earth, and the spectral radiance of a blackbody at 
6000° K. dose to the tentperature of the sun. On estimating the sky radiance one can 
assume that the radiation from the sun is diffusely scattered and hence the sky appears 
to have a uniform radiant energy. This kind of maximum brightness, which has an 
average radiance about 2 xlC' times the radiance of the sun. can be approximated by 
a bright cloud shown as the lower dashed curve in the figure. But normally the sky 
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does not scatter this large amount of radiation; a nominal value is about one-tenth of 
this maximum spectral radiance, shown as the solid curve on the left, named "clear 
sky". So we can describe the combined eflect of the thermal radiation emitted at the 
temperature of the earth and the scattered sunlight as combination of the two solid 
curves shown in this figure. 
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2. Atmospheric Emission 

a. General 

The atmospheric emissivity in the infrared depends on the amount of water 
vapor, carbon dioxide and other gases which give rise to absorption bands due to their 
vibrational and rotational transitions. Since the absorptivity equals the emissivity 
( Kirchhofi's law ). we can predict the atmospheric emission from the absorption 
spectra of these gases. 

b. Spectral Radiance of Clear Sky 

Sky temperature and elevation angle are the main factors determining the 
radiance of a clear sky. The elevation angle determines the length of the atmospheric 
path and hence the emissivity, w'hile the temperature determines the blackbody 
emission characteristics. 

The effect of the elevation angle on spectral radiance is presented in Fig 
3.12 [Ref. S]. At low angles of elevation ( near the horizon ) we can see nearly 
blackbody radiation at the ambient temperature. But as the elevation increases above 
the horizon and the atmospheric path length decreases the optical density decreases 
and the emissivity is no longer blackbody ( £ = 1 ) eveiywhere. In the regions 
centered about 6.3 pm { water vapor absorption band ) and about 15 pm ( carbon 
dioxide absorption band ) emissivity is very high, and the sky remains blackbody. even 
at the zenith ( 90°). But between these regions, in the region from about 8 pm to 13 
pm the emissivity is poor and the sky is no longer black. The general background 
emission in this region is due to the wings of the water vapor and carbon dioxide 
bands. The peak that we see at about 9.6 pm is due to ozone emission, which 
increases as the elevation angle decreases from the zenith. 

Fig 3.12, taken from [Ref S], shows spectra of clear sky radiance measured 
in September, at night, from a height of 11000 feet above sea level, at the ambient 
temperature of 8°C. The angles of elevation in sequence are 0°. 1.8°, 3.6°. 7.2°. 
14.2°, 30° and 90°. 

The effect of the ambient temperature is very strong for this kind of 
radiation, as shown in Fig 3.13 . This figure is also taken from [Ref. 8] and shows the 
zenith sky spectral radiance for two different ambient temperatures. The dashed curves 
are the spectral blackbody radiances corresponding to the ambient temperatures for the 
two separate measurements at 27.5 °C and 2.5 °C. These measurements were taken with 
different amounts of air mass, which affected the emissivity and hence the radiance of 
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Figure 3.12 Sky radiance for several elevation angles above horizon. 

the spectrum in the regions with poor einittance [Ref. S]. The low temperature 
measurement was taken at 14110 feet above sea level and the higher temperature 
measurement at 6000 feet above sea level. The negative values shown are due to 
measurement error. 

c. Spectral Radiance of an Overcast Sky 

Overcast skies have been found to be ver>' good blackbody radiators ( £ = 

1 ), having a temperature at low altitudes within one or two °C from the ambient 
surface temperature [Ref S]. 

As mentioned before, water in the infrared is a good absorber (see Fig 3.3) 
with a small rellectivity. Clouds mostly at lower altitudes consist of water droplets with 
radii from 1 j.im to about 10 |im at concentrations above 300 cm*'’ with temperatures 
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MICRONS wavelength 

Figure 3.13 Variation of zenitlt sky radiance with ambient temperature. 

not far from the surface temperature [Ref. 9]. So infrared radiation incident on a cloud 
water droplet is highly absorbed and partially reflected and transmitted. However due 
to multiple scattering from all neighbouring droplets in the cloud the transmittance of 
the cloud goes essentially to zero. 

Hence we can say that clouds in the infrared region are good absorbers and 
consequently good emitters, causing overcast skies to behave as good blackbody 
radiators. 
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d. Spectral Radiance of Clouds 

Clouds with thickness greater than 100 m and clouds with water droplet 
( or ice cr\’stal for higher altitudes ) counts greater than 0.5 g;cm~ do not transmit. 
Thus absorptivity, and consequently emissivity, is high, approaching one. which means 
that clouds show a ver>’ good blackbody behavior [Ref. 10]. 

Fig 3.14 [Ref 10] shows the emissivity of cloud layers versus wavelength 
and thickness AZ . .Also in Fig 3.15 [Ref 10] is shown how the integral emissivity (E). 
transmittance (T) and absorptance (.A) changes with thickness AZ. Here we can clearly 
notice that at thickness AZ above 60 m the transmittance goes to zero and above 100 
m the eniissivity goes to zero. 

However, looking at the spectrum of the under side of a cumulus cloud in a 
clear sky ( Fig 3.16 ) at an elevation angle 14.5° [Ref 8] we do not observe a 

blackbody spectrum . Since the emissivity of the cloud is expected to be high and since 

the sky is partially transparent in the region of 8 - 13^m for the elevation angle of 
14.5°, as shown before, we might expect to see clearly the emission of the cloud in that 
spectral region. From the figure we see that the emission in this region is characteristic 
of a blackbody at -10°C which is reasonable for the temperature of the cloud ( lower 
dashed line in figure ). What is really happening in this figure is that we have a 
combination of the atmospheric emission at the 6 and 15 j.im wavelength regions from 
the path below the cloud together with the -10° blackbody emission from the cloud in 
the wavelength region at 8 - 13 ftm . 

The effects of thin cirrus clouds in a clear sky modified situation is 

presented in Fig 3.18. Here we can see that the emission in region 8-13 J.im has been 

raised in comparison with the clear sky situation ( see Fig 3.13 ) , and also the ozone 
emission is reduced since it comes from a region higher than the clouds and is 
attenuated passing through them. 

Fig 3.17, from [Ref 8], shows spectra measured in June at sea level with 
an ambient temperature about 28°C. The elevation angles are 0°. 1.8°, 3.6°, 7.2°. 
14.5°. 30° and 90°. 

3. Scattering of the Sky Radiation 

The position of the sun relative to the area of the observed sky and the length 
of the scattering path due to the elevation angle affect strongly the amount of scattered 
radiation. We have a large amount of scattering close to the sun direction as well as an 
increase in scattered radiance with a decreasing elevation angle. Decreasing the 
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Figure 3.14 Spectral emissivity of cloud layers of thickness AZ. 




Figure 3.15 Integral eniissivity (E). absorptivity (.A.). transmittance (T). 

elevation angle we have a longer atmospheric path from which scattering can be 
produced. But verv' close to the horizon we have a reduction in the scattered radiance 
mainly because of water vapor absorption which reduces the elTective observation path. 

This kind of radiation is not important for the work done in this thesis since it 
occurs at short wavelengths. 
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Figure 3.16 Radiance of the under side of a cumulus cloud. 




Figure 3.17 Spectral radiance of sky covered with cirrus clouds. 
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IV. ATMOSPHERIC TRANSMITTANCE 



A. ATMOSPHERIC PROPAGATION 

Thermal radiant flux from a target has to pass tlirough the atmosphere of the 
earth before it reaches an infrared sensor. During this passage, the thermal radiation of 
the target is attenuated because of selective absorption due to several atmospheric 
gases and scattering by small particles suspended in the atmosphere, such as aerosols, 
rain, snow, smoke, fog and haze. Both absorption and scattering depend on wavelength 
(>.). 

In general we call atmospheric extinction the process by which the radiant flux 
from a target is attenuated passing through the atmosphere. This can be expressed as; 

t^O.) - e-'O-'R . 

where is the atmospheric transmission, erf/.) is the extinction coefficient and R is 

the range. The extinction coefficient is the sum of the absorption coefficient a(/.) and 
the scattering coefficient yfX) 

(T(X) = a(}.) + Y(X) . 

Both absorption and scattering cefficients are sums of molecular and aerosol 
components [Ref 2]. In the infrared region the absorption process is of greater 
importance (particularly in the 8 - 14 |im ) than the scattering process which is 
important in the visible and near infrared region. 

Water vapor, carbon dioxide and ozone are the most significant absorbing gases, 
having broad absorbing bands centered at the specified wavelengths: water vapor ( 2.". 
3.2, 6.3 ^im ), carbon dioxide ( 2.7, 4.3, 15 J.im ). ozone ( 4.8. 3.2. 6.3 [im ). The bands 
of carbon dioxide at 2.7 and 15 jim and water vapor at 6.3 (tm limit the atmospheric 
transmission in the infrared to two atmospheric windows at 3.5 to 5 |Um and at 8 to 14 
Hm. 

Infrared systems are designed to operate through these windows, and it is very 
important for the performance of our system to have a prediction of the atmospheric 
attenuation for the wavelength and the environmental conditions in which we are 
operating. This prediction process is a compheated problem and has been faced by 
computer modeling codes. 



Such niodels are the high resolution models ( HITR.A.N) and the low resolution 
model ( LOWTR.AN). The high resolution model codes provide direct computation 
with high resolution single frequency results. They are used to predict atmospheric 
attenuation for laser radiation. The low resolution model codes are used for a broad 
band prediction where the thermal systems work. 

LOWTR.AN code was used to predict atmospheric transmittance for the work 
done in this thesis and is considered in the next sections of this chapter. 

B. LOWTRAN 6 CODE 

1. General Description 

LOWTR.AN 6, the latest edition of the LOWTfCAN simulation codes, is a 
FORTR..AN computer curve-fitting program, relating empirical precomputed 
attenuation data to simple meteorological and environmental inserted parameters. It 
was designed to calculate average atmospheric transmittance and thermal radiance 
over a specified broad band, ranging from 359 to 40000 cm'^ ( 0.25 to 28.57 jim ) with 
a resolution of 20 cm'^ averaging bandwidth in steps of 5 cm'^. 

The average transmittance is a function of the total extinction, the product of 
the transmittances due to water vapor line absorption, water vapor continuum 
absorption, uniformly mixed gases line absorption, nitrogen continuum absorption, 
aerosol absorption, aerosol scattering and molecular scattering. A consequence of the 
averaging over a finite spectral band is that the rapid fluctuations of the atmospheric 
transmittance in that band vanish. 

Generally we can say that the LOWTR-AN code is a fast but not ver>' 
accurate computer simulation code. 

2. LOWTRAN 6 Input Data 

The LOWTR.AN 6 code is activated by submitting four cards of sequential 
data. Each card provides data to be selected and asks for data values to be inserted as 
follows: 

The first card provides a selection of .six atmospheric models. These models 
are the tropical atmosphere, midlatitude sunuiier, midlatitude winter, subarctic summer, 
subarctic winter and the 1962 U.S. standard atmosphere. Also in this card we have to 
indicate the type of atmospheric path ( horizontal, vertical or slant ), and determine the 
execution mode ( transmittance - radiance ). 



36 



The second card provides a selection of the type of extinction for nine 
bountaiy- layer aerosol models and gives a default meteorological range. These models 
are the rural (vis = 23 km), rural (vis = 5 km), naw maritime (computes visibility), 
maritime (vis = 23 km), urban (vis = 5 km), tropospheric (vis = 50 km), advection 
fog (vis = 0.2 km), radiation fog (vis = 0.5 km) and a user defined model (vis = 23 
km). Also here we have to specify the air mass character for the navv' maritime model 
and some meteorological conditions. 

The third card defines the geometrical path and we have to insert some 
parameters like initial and final altitude, initial zenith angle and length of our path. 

In the fourth card we have to specify the spectral range over which we want 
data for atmospheric transmittance or radiance and the frequency increment over 
which we want our data to be averaged. 

•A. detailed description of the functions of these cards and the required input 
data is given in [Ref. 11]. 

C. LOWTRAN 6 CODE USED IN THIS \VORK 

For the work done in this thesis the LOWTRAN 6 code was used to find the 
atmospheric transmittance for the different meteorological conditions we experienced . 

For our experimental conditions, mostly overcast skies and ambient temperatures 
about 10°C to I6®C we selected the midlatitude winter ( 45° N. Januarx' ) atmospheric 
model for horizontal path. Also for the purpose of this work the na\w maritime 
bountaiy layer aerosol model was selected. .And thus because this model takes into 
account the existence of three components in determining the aerosol population. First 
is the coniinenial component which is sensitive to the amount of time the air mass has 
spent over the sea; Second is the stationary component which depends on the wind 
history and represents aerosols formed by high winds and whitecap conditions, and the 
third one is the fresh component which is characterised by the current wind speed and 
consists of whitecap produced droplets. For our geometrical path we inserted target 
altitude 10 m and horizontal path length 650 m. Our spectral range was 710 - 1250 
cm'^ (8-14 ^m ) with frequency increment 5 cm'^ 

The results of these calculations are discussed in Chapter VH. 
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V. FAMILIARIZATION WITH THE EQUIPMENT 



A. GENERAL 

As mentioned in the introduction, all our data for this work were collected, 
stored and processed using the AGA Thermovision 780 s}'stem. supplemented with a 
digital microcomputer. 

The system is designed for both thermal Imaging and image analysis. It 
combines real infrared scanning with thermal image measurement capabilities and is 
supplemented by a microcomputer for digital image recording, storing and processing. 

The particular system we used was a dual scanning system operating in the 3 - 
5.6 fim and -8 - 14 fiin windows with two different detectors. For the purpose of this 
experimental work only the 8 - 14 fim scanner was used. 

B. AGA THERMOVISION 780 

1. Basic Description 

The AGA Thermovision 780 system consists of a real time scanner, a black 
and white monitor chassis display and a digital microcomputer. 

The real time infrared scanner is the unit that converts the thermal radiant 
energy emitted from a target to electronic video signals. These signals are amplified 
and transmitted through a cable to the display monitor chassis where with further 
amplification they produce a thermal image on a black and white screen. The amplified 
signal can also be fed to a color monitor through an interconnecting device in order to 
have a real time thermal image presentation. Finally video and data signals are 
transmitted through an interconnecting link to a microcomputer. Here the thermal 
image is desplayed in color and can be recorded, stored on a data disk and processed 
for temperature evaluation using color codes. 

2. The Infrared Scanner 

There are two kinds of scanning unit used in this system, the short wave (SW) 
scanner that operates in the 3 to 5.6 fim spectral band and the long wave (L\V) scanner 
that operates in the 8 to 14 jim band. Scanner systems are designed with single or dual 
scanner units which cover both bands simultaneously. Our scanner as we have already 
mentioned was a dual type. Such a scanner system is shown in Fig 5.1. 

Each scanner unit consists of the following parts : 
a. Electro-optical scanning mechanism 
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b. Infrared detector 

c. Liquid nitrogen dewar for cooling the detector 

d. Control electronics and preamplifier 

A simplified block diagram of the scanner unit is shown in Fig 5.2. 




Figure 5.1 Dual scanner system. 




Figure 5.2 Simplified scanner unit block diagram. 
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Electromagnetic thermal energy emitted from a target is focussed by infrared 
lenses into a vertical prisnt. The lenses for the short wavelength band (SW) are made of 
silicon and those of the long wavelength band (LW) are made of germanium. Both 
kinds are covered with proper antireflecting coatings. The vertical prism is rotated at 
ISO rpm and its optical output is passed through a horizontal prism which rotates at 
ISOOO rpm. 

Rotating slotted discs control the rotation of both prisms and are 
electronically connected to the horizontal and vertical triggering circuits to provide 
horizontal and vertical triggering pulses for the monitor. The horizontal and vertical 
prisms are synchronized in such a way that four fields each having 100 horizontal 
scanning lines produce one interlaced frame. Of these lines only 70 lines per field or 
280 per frame are used as active imaging lines. 

The output from the horizontal prism passes through a set of relay optics 
containing a selectable aperture turret with eight aperture discs between f l.S and f 20 
(f l.S disc was selected for our case ). It then passes through a filter unit ( no filter for 
our case ) and finally is focussed on a single detector located on the wall of the dewar 
chamber. The detectors are a photovoltaic Indium-Antimonide (InSb) for the short 
wavelength band and a photoconductive Mercur\’ - Cadmium - Telluride (HgCdTe) for 
the long wavelength band. Both detectors are cooled with liquid Nitrogen which 
maintains the chamber at a temperature of -196° C. 

The signal output from the detector is preamplified within the scanner unit 
and the video signal produced is applied through a cable to the black and white 
monitor chassis [Ref 12]. 

The arrangement of the electro-optical components of a scanner unit is shown 
in Fig 5.3 [Ref 12]. 

3. Black and White Monitor Chassis 

The video output signal from the scanner unit is applied to the input 
connector of the black and white chassis. This signal after a further amplification and 
processing within the monitor is applied to the display screen. 

Outputs from this monitor chassis are used for a color monitor connection 
and the interconnecting link for the microcomputer connection with the system. 

A Black and White monitor chassis is shown in Fig 5.4 [Ref 12]. 
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Figure 5.3 Arrangement of electro-optical components in a scanner. 




Figure 5.4 Black and White monitor chassis. 
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4. Technical Data 

Technical data for the AGA Thermovision 780 were found in [Ref. 13] and are 
presented in Table 1. 

5. Digital Image Processing System 

This system supplements the AGA Thermovision 780 system and consists of 
an interconnecting link and the microcomputer which utilizes the systems software 
DISCO 3.0 

a. Interconnecting Link 

The interconnecting link is a digital on-line interface for direct interfacing 
between the Thermovision 780 and the IF 800 model 20 microcomputer. The main 
purpose is to digitize the image, hold the image data during the time required for 
transfer from the Thermovision to the microcomputer, and control the data 
transferring speed. 

b. IF SOO Microcomputer 

IF 800 is a business level personal computer which uses a Z-80A . 4 MHz 
clock microprocessor for its CPU. It comprises a keyboard, a printer, a display unit 
and a mini floppy disk drive. 

6. DISCO 3.0 Software Description 

The software used by this system consists of three main programs the DISCO, 
the IMAGE PROCESSING and the UTILITY programs. Various subprograms as 
shown in the DISCO 3.0 block diagram Fig 5.5 [Ref. 14] supplement the main 
programs. 

The main program DISCO was the one mostly used for this work to record, 
store and process our data. Using this program infrared images were transferred from 
the Thermovision 780 system to the microcomputer. The subprogram FI Record 
Image was used to record these thermal images and display them in color on the 
computer screen. Also this subprogram was used to enter or alter our parameters 
(emissivity, atmospheric transmittance etc.), to select and enter the scanner 
identification characteristics used and store the images on a data disk. The subprogram 
F2 Evaluate Image was used to evaluate the tentperature measurentent of the image 
using color scaling with tables of 3 to 36 different color combinations. Each color was 
assigned a particular temperature range. The standard table of the system uses 8 colors 
and was the one used in our case. This subprogram was also used to display 
histograms with a frequency distribution of the temperature values, to display single 
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TABLE 1 

AGA THERMOVISION 7S0 TECHNICAL CHAR.ACTERISTICS 

FREFORMAHCE 



Spectral Range • 


3-5.6 |tm and/or 8-14 jim 


Frame Rate • 


6.25 sec' ^ 


Field Rate • 


25 sec”^ 


Interlace i 

REPLACEABLE FORE OPTICS 


4 j 1 


FOV Azimuth * i 


7* 


FOV Elevation * j 


7«* 


IFOV Azimuth : 


1. 1 mrad 


IFOV Elevation : 


1. 1 mrad 


MEAT ! 


0. 12»c at 22«c 


HDT t 


< 0. l®c 


Dynamic Range : 

OPTICAL DATA 


-20 to 900«c 



Effective Aperture Area : 24 

Aperture Diameter : 5.5 cm 

Effective Focal Length ! 9.9 cm 



f /number 
DETECTOR 


: 1. 8 


Type SW j 


Photovoltaic InSb 


LW j 


Fhotoconduct i ve HgCdTe 


Number of Elements : 


1 


Peak Wavelength SW • 


5 jim 


LW ! 

COOLING SYSTEM 


10 |tm 


Lquid Nitrogen : 


-196«c 



* Note! The system can 


take five different fore lenses with 



FOV 3. 5X3. 5, 7X7, 12X12, 20X20. and 40X40. The l*»n3 
we us»»d was the one with FOV 7X7. 



43 



DISCO 3.0 
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Figure 5.5 Block diagram of DISCO 3.0 software. 



isotherms, magnify the image, transfer a new image on the display from those which 
had been stored on the data disk and print data such as temperature profiles using the 
system printer. For the last case our system had some limitations since the system's 
maximum capabilities are available only on a color printer which was not available. 

C THERMAL MEASUREMENT TECHNIQUES 

1. Introduction 

As we know the radiant flux density received by an infrared system does not 
have a linear relationship to the object temperature. It is always mixed with rellected or 
emitted background radiation and is attenuated by the atmosphere. Therefore, certain 
calibration and correction procedures have to be applied to a particular measurement 
in order to have an accurate temperature evaluation. 

The AG.A Thermovision 780 system measures thermal radiation between 
spectral limits as mentioned before. The numerical value of the received rad'iation is 
called the THERMAL VALUE and as measured in Isotherm Units (lU) which is an 
arbitrary unit system. This thermal value is proportional to the received photon flux 
density, but is not proportional to the temperature of the target. Therefore, we need a 
calibration function as a graphical curve or a computer program to relate the thermal 
value ( isotherm units ) with the temperature of the target. 

2. Calibration Curves 

Calibration curves for all aperture sizes used by the system in the long 
wavelength ( LW ) spectral band, in which we worked, are shown in Figs 5.6 through 
5.8 [Ref 12j. These curves have been obtained assuming emissivity C = 1 (blackbody 
radiator ) and neglecting atmospheric damping = 1. For situations dilferent from 
the above it is obvious that we have to apply correction techniques. 

These calibration curves are accurately described by the following equation 
taken from [Ref 12] 

A 

I = , (eqn5.1) 

C e^^ - 1 

^ I 

where I is the thermal value in isotherm units ( lU ) for the absolute temperature (T) 
in °K. A,B,C are calibration constants that depend on the aperture, filter, scanner etc. 

This equation was used by the microcomputer program to calculate the 
temperature of the target for a known thermal value after the calibration constants 
•A.B.C for the LW shown in Table 2 had been inserted. 
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>20 -10 0 10 20 30 40 50 

temperature TC 



Scanner ? LVD9 4011 Date: 62"*09“'27 Operator? 

Detector ? C 1739 Ambient temp ? *C 

Filter ; Rel. humidity ? ^ % 

Lens^^^ ^ 3105 Object distance ? ^ m 



Figure 5.6 AGA LW calibration curves from -2QoC to 50oC. 
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'Scanner : LVDB 4011 Date : 82-09-27 Operator : 6. A 

Detector: G 1739 Ambient temp: 24 'c 
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7- 3105 
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m 



Figure 5.7 AGA LW calibration curv^es from -20 °C to 150°C. 
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Figure 5.8 AG A LW calibration curves from 0<>C to lOOOoC. 



TABLE 2 



AGA THERMOVISION 780 
LW CALIBRATION CURVE CONSTANTS 



INDIVIDUAL CALIBR.ATION OF 7S0 DUAL BBAR 



SERIAL NUMBERS 

SCANNER :LWDR4011 
DETECTOR ; G 1739 
FILTER : NOF 

LENS :7 3105 

CALIBIUATION CONDITIONS 

AMBIENT TEMPEIUATURE : 24 C 

RELATIVE HUMIDITY : 55% 

OBJECT DISTANCE : 1.0 m 

CALIBIUATION DATE ; 82 - 09 - 27 



CAL1BR.ATION CURVE CONSTANTS 



APERTURE 


A 


B 


C 


1.8 


- 3581 


1506.49 


- 0.436 


2.5 


- 4060 


1569.30 


- 0.759 


3.6 


- 6514 


1629.70 


- 1.824 


5.1 


5420 


1610.70 


2.796 


7.2 


1123 


1606.54 


1.098 


10.0 


584 


1604.29 


0.885 


14.0 


306 


1610.60 


0.773 


20.0 


195 


1650.51 


0.767 
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3. Basic Thermal Measurements 

There are two basic methods we can use to measure temperature using the 
thermovision 780. The direct, and the relative measurement method. 

In the direct method of measurement the system built in capabilities (clamping 
and temperature compensation) permit measurement of the temperature of an object 
without using any other reference source of radiation. 

In the relative measurement method the radiation of an external reference 
source with known temperature and entissivity is compared with the radiation received 
from the object to obtain the object temperature. 

The direct method is easier to use but not as accurate as the relative method, 
which, however, has the disadvantage that it needs an accurate reference source to be 
available. In our case the direct method was used since we didn't have any reference 
source with known temperature and emissivity located near the target. 

4. Direct Measurement Procedure 

In this subsection we will describe briefly the procedure for the direct 
measurement method, using the Thermovision first alone and then supplemented with 
a microcomputer as was done to collect data for this work. 

To aid in understanding this measurement procedure, the monitor chassis 
front panel controls and indicators are presented in Fig 5.9 [Ref. 12]. Subsequently a 
short description of the operation and function of the controls and indicators referred 
to in this measurement procedure is given. 

a. THER.M.A.L RANGE (item 9 ); This is a nine position switch calibrated 
between 2 and lUUU isotherm units. This switch selects the thermal span of 
interest. 

b. THERMAL LEVEL (item 10 ): This control sets the thermal level of the image. 

c. PICTL'RE MODE "NORM.VL" ( item 12b ): It selects the normal imaee 
displav, where the cooler object temperature appears as black and the hotter as 
white.' 

d. ISOTHERM LEVEL 1 (item 13 ); This control turns on the isotherm function 
and identifies the discrete thermal levels on the target. 1 he isotherm function is 
indicated by saturated white areas, which represent areas of the same 
temperature) .A, marker on the vertical isotherm scale indicates the relative 
position of the selected level. 

e. ISOTHERM SCALE (item 16 ): This is a scale calibrated from -0.5 to +0.5. A 
reading of the isotherm marker multiplied bv the thermal range gives the 
relative thermal value in isotherm units. 

For both cases either with the thermovision itself or supplemented by the 
microcomputer, starting with picture mode in NORMAL, we have to adjust the 
THERMAL ICANGE and LEVEL controls to obtain a good thermal picture. The 
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Figure 5.9 Monitor chassis front panel controls and indicators. 



picture was shown on the Black'White monitor and on the color monitor used. The 
measurement procedure is shown schematically in Fig 5.10 [Ref. 12]. First with the 
THERMAL R.AXGE control switch we select the thermal span of interest in isotherm 
units ( 2 to 1000 ). about the thermal level of the image. In our case thermal ranges of 
2. 5 and 10 isotherm units were selected. Then with the THER.MAL LEVEL control 
knob we set the thermal level of the image. .At the point at which a satisfactory 
thermal picture is obtained the corresponding thermal value "L" of the image in 
isotherm units, is measured (i.e., the thermal value of the isotherm scale at zero point. 
Fig 5.10). 

Using the thermovision alone, the control ISOTHERM LEVEL I has to be 
adjusted to brighten up the point of interest on the target in view. A marker on the 
vertical scale indicates the relative position of the selected level. The reading of the 
isotherm marker multiplied by the thermal range gives the relative thermal value "i" in 
isotherm units, shown in Fig 5.10. .Adding the two values "L" and "i" we have the 
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NOTE; 




thermal LEVEL = L 
RELATIVE THERMAL VALUE = i 
measured THERMAL VALUE = I^ =L + i 



(°C) 

8LACKB0DY TEMPERATURE 



TEMP Range 



Figure 5.10 Direct measurement method. 

measured thermal value 1'^ for the particular point of interest. Under the assumed 
simplified conditions, having = 1 and e = 1 , the calibration curve can be used, as 
in Fig 5.10, to translate the measured thermal value 1'^ to temperature. 

Using the microcomputer as was done to collect data for this work there was 
no need to go through the isotherm procedure. The computer does not use isotherms 
but discrete thermal values. Therefore, we had the isotherm marker at zero. The 
subprogram FI Record Image was used to record not only the thermal image picture 
but also the corresponding themial value and thermal range. These values were used by 
the subprogram F2 Evaluate Image to evaluate and present on the display the 
temperature distribution of the recorded image in color scaling, using the standard 
eight color combination table of the system. 





D. COMPLEX THERMAL MEASUREMENTS 

1. Introduction 

As mentioned before, the AGA Thermovision 780 system considers an ideal 
situation for the values it measures. But as we know there are several factors that 
influence these measurements. So the true object measurement must be derived bv a 
calculation in which we must take into account all the factors that influence the 
received thermal radiation. 

2. Variables Affecting the Measurement 

The most important factors that aflect the measured thermal radiation are the 
following: 

a. Emissivity of the target 

b. Background surrounding the target 

c. Target opacity 

d. Target size 

e. Atmospheric effects 

Most of the objects in our environment are not blackbodies and hence their 
emissivities are less than one. For our case the emissivity of the R V "POINT SLR" 
was found to be about 0.95 as will be explained later. The influence of the background 
was discussed in detail in Chapters II and III and the atmospheric eflects were 
discussed in Chapter IV. 

The target opacity has to be taken into account because some targets are 
more or less transparent to infrared radiation ( glass, plastics, clouds, etc. ). Radiation 
received from objects like these includes radiation transmitted from the background. 
Also the emissivity of the target is aflected and is hard to estimate. For our case the 
target was considered as an opaque Lambertian surface. 

The size of the target is another factor to be taken into account. Since the 
detector subtends a specific solid angle, if the target is not big enough to cover this 
angle the detector receives radiation from the target background. This will tend to 
aflect the temperature difference indicated between the object and the background. In 
our Case the target was not big enough to cover the field of view of a 7°X7° lens at a 
distance 650 m. For most measurements it occupied 13 to 2 5 of the display screen. 

3. Exact Measurement Formula Derivation 

The thermal radiation that reaches the detector of our system is the sum of 
the target-emitted thermal radiation, the radiation reflected by the target and the 
radiation emitted from the atmosphere as shown in Fig 5. 11. 
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Figure 5.1 1 Radiation conditions in complex measurement. 

If we assume an opaque object at temperature (T^) and emissivity ( ), this 

object emits photon flux CqS^. is the photon flux of a blackbody at the same ( T^) 
temperature. Since the above llux is attenuated by atmospheric transmittance ( ) the 

object radiation reaching our detector is 

.Also assuming that all surrounding surfaces are of the same ambient 
temperature ( ) which determines that the emissivity is close to one ( = 1 ) we 

can say that the radiation reflected by the object is PqS,^. Sy is the photon llux of the 
surroundings emitting as a blackbody ( Cy = 1 ) at temperature ( Ty ) and is the 
reflectance of the object. From the first assumption that the object is opaque we get 
Pq = * ■ ^0 since the above radiation is attenuated by the atmosphere the 

reflected radiation reaching our detector is "y( 1 - )Sy. 

Finally the atmospheric radiation is given by ( 1 - Ty where is 

the photon flux due to atmospheric emission and scattering. Hence we can write that 
the photon flux S'^ reaching the detector is 

S o - Cyi-ySy + Ty(l - c^jSy + ( 1 ' ^ (j( 1 S ^ (cqn 5.2) . 

This photon llux relation can be converted into a thermal units relation since, 
as we mentioned before, the thermal value ( lU ) measured by the .AGA Thermovision 
7S0 is proportional to the photon flux reaching the detector. Therefore we can write 
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the relation I — C X- S where I is the thermal value. C the proportionality constant 
and S the photon flux. 

Substituting S = I,C in Equation 5.1 we have the same equation in thermal 

values 

^ o = ^a^o^o + - ^o)I« + - V latm 5.3) . 

This equation can be solved for to find the exact temperature through the 
calibration curves or can be solved to find the emissivity if we know the exact 
temperature of the object and consequently from the calibration curves. Therefore 
from Equation 5.3 we have 



I 



lo = 



0 I 11 

--(--l)^a--(--niatm 



^0^0 >^o 



^0 



(eqn 5.4) or 



Gq = (eqn 5.6) . 

^a^^o * 

For our case the computer program calculated for us the exact temperature 
using Equation 5.3 after we had inserted certain values for the parameters emissivity. 
atmospheric transmittance, ambient temperature and temperature of the atmosphere. 
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VI. EMPIRICAL CALIBRATION OF AGA THERMOVISION 780 



A. CALIBRATION PROBLEM 

Trying to process the thermal images we had already recorded, we realized that 
the temperatures the computer was calculating were much higher than expected for the 
particular scene. 

It was obvious that the Thermovision 780 was measuring thermal values ( lU ) 
higher than the values the calibration curves were giving for certain target 
temperatures. 

To solve the problem we had to derive new calibration curves to use, for the 
exact temperature evaluation of the stored thermal images. At that time we had 
difficulty deriving new calibration constants to be inserted in the computer program. 
Therefore, we tried to find a relation between the exact temperature of a blackbody 
and the temperature the computer program was calculating using the faulty calibration 
curves for the blackbody thermal values measurements. In this way we had an 
empirical relation as a curve function, accurate enough to correct the temperatures the 
computer program was calculating using the faulty calibration curve on processing the 
thermal images we had stored. 

B. THE CALIBRATION PROCESS 

1. General 

For the calibration procedure we required a blackbody working at about the 
temperatures we had when the images were recorded. The air temperature at that 
period was 9.7 - 11.3°C, the sea surface temperature was 11.2 - 12.8°C and the 
laboratory ambient temperature, where the calibration procedure was carried out, was 
always 18 - 19® C. Therefore we needed to have a blackbody source to work below the 
ambient laboratory temperature in order to get the calibration relation wc wanted. 

2, Blackbody Source 

For the purpose of this work we built our own blackbody source to use in our 
calibration procedure. As we have already mentioned we wanted this source to work 
at about the expected target temperatures, which vv'ere below" the laboratory ambient 
temperature. 
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The core of this blackbody source was a cylinder 2S cm long and 7.5 cm wide 
made of thick ( 0.5 cm ) aluminum. One end of this cylinder was cut at an angle ab'out 
60® from horizontal and it was covered with a plate of the same thickness of aluminum 
for minimizing rellected radiation escaping out of the cavity. The other end was cut at 
90° and it was covered with a same thickness aluminum plate which had a hole of 2 
cm diameter. The interior was painted with mat high emissivity black paint. Using a 
nichrome \wire wound uniforntly around the cylinder we could heat our source and by 
pouring liquid nitrogen inside the cylinder we could cool it below ambient temperature. 
The whole source was covered with a kind of fiber glass insulating material and four 
thermocouples were used to check the temperature stability of the core cylinder and the 
edge plates. In Fig 6.1 a schematic representation of this blackbody source is shown. 

This source was well tested showing ver>' good blackbody behavior, stabilizing 
and maintaining the temperature we needed quite fast. 




Figure 6.1 Blackbody source used for calibration procedure. 

3. Calibration 

We wished to find a relation between the exact temperature of an object and 
the temperature calculated by the computer program using the wrong calibration 
curx’es for the isothermal values measured by the Thermovision 7S0. 
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For this reason our blackbody was set to certain temperatures about the 
expected target temperature at the time of the experiment. After these temperatures 
were stabilized a measurement with the Thermovision 7S0 was taken from a distance of 
about 2 m. The thermal value was transferred to the computer and the thermal image 
was recorded, to be processed for temperature evaluation. For the distance of 2 m the 
atmospheric attenuation was neglected ( ~ 1 .)• The above procedure was repeated 

for all three thermal ranges ( 2. 5, 10 ) used to collect data. 

The above calibration temperature measurements for the three different 
thermal ranges are shown in the following Tables 3,4.5. Looking at these 
measurements we can notice that there is a linear relation between the blackbody exact 
temperature and the temperature the computer was calculating. Using the method of a 
least squares linear regression we found the equation of the best fitting straight line. 
This line was the calibration curve function we wanted to correct the temperatures of 
the stored thermal images. 

Plotted calibration data measurements and the best fitting straight line for all 
three thermal ranges are shown in the following Figures 6.2 to 6.4. 

C. CONCLUSIONS 

From the calibration curves given in the previous chapter we can see that the 
portion of the curves in the temperature region used for the calibration procedure is a 
straight line. The linearity of that region was the confirmation that the linear relation 
we were getting was reasonable. 

The corrections we were getting from these straight lines were accurate to within 
an error of about 5%. The thermal values measured by the thermovision were 
accurate to within half a thermal unit. This caused an error of 1.5 % in the calculated 
temperatures. The variance of the best fitting staight line for each thermal range caused 
a 3 % error in the thermal range 2, 3.5 % in the thermal range 5 and 3.7 % in the 
thermal range 10. Therefore, the total correction error was within 4.5 % in the thermal 
range 2, 5 % in the thermal range 5 and 5.2 % in the thermal range 10. Of course we 
have to keep in mind that these calibation functions are only valid for the temperature 
region about our target temperatures ( 9 to IS^C ). 
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TABLE 3 



THERMAL R.ANGE 2 

CALIBR.ATION TEMPER.ATURE MEASUREMENTS 



TEMPER.ATURES IN DEGREES C 



BLACKBODY AGA 

EXACT COMPUTED 


DIFFERENCE 


3.0 


16.1 


13.1 


4.2 


16.5 


12.3 


5.0 


17.2 


12.2 


5.6 


17.8 


12.2 


6.3 


18.2 


11.9 


7.0 


18.7 


11.7 


8.0 


19.7 


11.7 


8.5 


20.1 


11.6 


9.0 


20.5 


11.3 


9.5 


21.2 


11.6 


10.0 


21.8 


11.8 


11.0 


22.5 


11.5 


12.0 


23.4 


11.5 


12.5 


23.8 


11.3 


13.0 


24.5 


11.5 


13.5 


25.1 


11.6 


14.0 


25.7 


11.7 


14.5 


26.1 


11.6 



TEMPER.ATURES IN DEGREES C 

BLACKBODY AGA DIPTERENCE 
EXACT COMPUTED 



15.5 


26.9 


11.4 


16.0 


27.3 . 


11.3 


16.5 


27.6 


11.1 


17.0 


28.0 


11.0 


17.5 


28.3 


10.8 


18.0 


28.9 


10.9 


18.5 


29.4 


10.9 


19.0 


29.8 


10.8 


19.5 


30.2 


10.7 


20.0 


30.7 


10.7 


20.5 


31.1 


10.6 


21.0 


31.6 


10.6 


21.5 


32.0 


10.5 


22.0 


32.4 


10.4 


22.5 


33.0 


10.5 


23.0 


33.6 


10.6 


23.5 


33.9 


10.4 



MEAN DIFFERENCE : 11.03 
STANDARD DEVIATION ; 1.724 
VARIANCE : 2.SS7 

BEST FIT STR.AIGHT LINE : Y = (1.0238) X - 11.62 
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Figure 6.2 Tlicrmal Range 2 calibration temperature measurements. 



TABLE 4 



THERMAL R.ANGE 5 

CALIBRATION TEMPEFLATURE MEASUREMENTS 



TEMPEIUATURES IN DEGREES C 



.ACKBODY 


AGA 


DIFFERENCE 


EXACT 


COMPUTED 




0.0 


14.1 


14.1 


1.0 


15.1 


14.1 


2.0 


15.8 


13.8 


3.0 


16.8 


13.8 


4.0 


17.6 


13.6 


5.0 


18.2 


13.2 


6.0 


18.9 


12.9 


7.0 


19.2 


12.2 


8.0 


20.4 


12.4 


9.0 


21.7 


12.7 


10.0 


22.3 


12.3 


11.0 


23.0 


12.0 


12-0 


24.2 


12.2 


13.0 


25.0 


12.0 


14.0 


25.7 


11.7 



TEMPER.ATURES IN DEGREES C 

BLACKBODY AGA DIFFERENCE 
EXACT COMPUTED 

16.0 26.7 10.7 

17.0 27.6 10.6 

15.0 2S.2 10.2 

19.0 29.6 ■ 10.6 

20.0 30.5 10.5 

21.0 31.0 10.0 

22.0 32.3 10.3 

23.0 32.6 9.6 

24.0 34.2 10.2 

25.0 34.8 9.S 

26.0 36.5 10.5 

27.0 36.9 9.9 

28.0 37.6 . 9.6 

29.0 38.5 9.5 



ME.AN DIFFERENCE : 11.55 
STANDARD DEVIATION ; 1.531 
VARI.ANCE : 2.265 

BEST FIT STR.AIGHT LINE : Y = (1.1956) X - 16.646 
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1 igure 6.3 1 hermal Range 5 calibration temperature measurements. 



TABLE 5 



THERMAL RANGE 10 

CALIBR.AT10N TEMPER^ATURE MEASL'REMETNS 



TEMPER.ATURES IN DEGREES C TEMPER.ATURES IN DEGREES C 



BLACKBODY AGA DIFFERENCE BLACKBODY AGA DIFFERENCE 



EXACT 


COMPUTED 


0.5 


14.2 


1.5 


15.0 


2.5 


15.6 


3.5 


16.8 


4.5 


17.7 


5.5 


17.9 


6.5 


18.6 


7.5 


19.4 


8.5 


20.1 


9.5 


20.6 


10.5 


21.0 


11.5 


22.9 


12.5 


24.1 


13.5 


25.2 


14.5 


25.8 





EXACT 


13.7 


16.5 


13.5 


17.5 


13.1 


19.5 


13.3 


20.5 


13.2 


21.5 


12.4 


22.5 


12.1 


23.5 


11.9 


24.5 


11.6 


25.5 


11.1 


26.5 


10.5 


27.5 


11.4 


28.5 


11.6 


29.5 


11.7 




11.3 





COMPUTED 




27.2 


10.7 


27.5 


10.0 


29.5 


10.0 


30.6 


10.1 


31.2 


9.7 


32.4 


9.9 


32.8 


9.3 


34.3 


9.8 


35.0 


9.5 


35.6 


9.1 


36.9 


9.4 


37.5 


9.0 


38.6 


9.1 



MEAN DIFFERENCE : 11.00 
STANDARD DEVIATION : 1.4S9 
VARIANCE : 2.141 

BEST FIT STR.AIGHT LINE : Y = (1.1807) X - 15.656 
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I’igiire 6.4 1 hermal Range 10 calibration temperature measuieinents. 



VII. SIGNATURE MEASUREMENTS 



A. GENERAL 

All signature measurements were taken during mid May 1986, on the following 
dates and time periods. 

1. .May 14 09:00 - 09:25 

2. .May 18 11:57-13:07 

3. .May 21 00:05 - 00:13 

4. .May 21 17:18-17:37 

5. .May 22 00:48 - 01:03 

Thermal imaging data were collected for the spectral range 8-14 pm using the 
LW detector of the AG.A. Thermovision 780. The thermovision system and its 
supporting computer were set up in a laboratory building at the eastern side of the 
HOPKINS .MARINE ST.VTIO.N grounds, at a window facing directly over the shore 
line toward the closest buoy. The laboratory is about 8 - 10 meters above sea level. 

The target was the oceanographic ship R/V POINT SUR which was on an 
operational oceanography student cruise. During those periods in which it was made 
available for signature measurements the ship maneuvered in the vicinity of the off 
shore ( 635 meters ) buoy, displaying various aspects as requested. The distance 
between the the ship and our detector was about 650 m. At that distance, the closest 
the ship could reach, using a 7<>X7® FOV lens we had an image covering almost 2 3 of 
the screen. Six thermocouples installed about the superstructure were used to have a 
calibration of the actual ship temperature. The thermocouples were sampled at 
10-seconds intervals and the output voltages were recorded, to be later converted into 
temperatures. 

The environmental conditions for the time period in which the measurements 
were taken are shown in the data presentation tables. Generally for these periods the 
area was under the influence of a high pressure system, with North West winds and 
coastal fog, especially during the morning hours. 

Thermal images for different ship aspects, were recorded and stored during time 
periods in which the ship was available for us and the environmental conditions 
allowed. All data were recorded assuming the ideal conditions appropriate to the 
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Thermovision system measurements, that the target was a blackbody source and 
atmospheric attenuation was neglected. 

Fig 7.1 is a map of the Monterey Peninsula showing the area where the 
signature measurements were taken. 




Figure 7.1 Area where the signature measurements were taken. 
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B. PARAMETERS FOR TEMPERATURE EVALUATION 

As mentioned in Chapter five the processing of the recorded and stored thermal 
images was done using the DISCO 3.0 computer program. For temperature evaluation 
the program asks for certain parameters to be provided. These parameters, as we have 
discussed before in the derivation of the exact measurement formula, in Chapter five, 
are the atmospheric (air) temperature, the ambient temperature, the atmospheric 
transmittance and the target emissivity. 

1. Atmospheric Transmittance 

The atmospheric transmittance was calculated using the LOWTICAN 6 code. 
For our experimental conditions, temperatures about 10 to looC, we selected the 
midlatitude winter (45°. January’ ) atmospheric model for a horizontal path. The Navy 
Maritime was selected for the boundary layer aerosol model, discussed in detail in 
chapter four. For the geometrical path we inserted target altitude 10 m and horizontal 
path length 630 m. The spectral range was set at S - 14 jim ( 710 1250 cm'^ ) with 

frequency increment 5 cm'^ 

The meteorological data needed, mainly wind speed and wind direction, for 
determining the air mass character, continental influence, were taken from 10-minute 
logged meteorological data recorded from instruments stationed on the ship. Such data 
weren't recorded for the last two observation periods on May 21 ( 17:18 - 17:37 ) and 
May 22 ( 00:48 - 01:03 ). For these periods meteorological conditions were 
approximated using observations taken from the meteorological station at POINT 
PINOS, the closest station in the area of the experiment ( .see Fig 7.1 ). Since the 
direction of the wind was from North West the continental influence as we can see 
from the map in Fig 7.1 was not very strong. So on a scale from 0 ( open ocean ) to 10 
( maximum continental influence ) we inserted mass charaeter scale from 2 to 3 . 

In Table 6 are shown the meteorological data used, the air mass character 
scale ( ICSTL ) and the calculated atmospheric tranmittance for each observation 
period. 

2. Emissivity 

The emissivity of the ship surface was found by evaluating the emissivity of 
two metallic pieces taken from different parts on the ship. These specimens were 
painted with slightly different white paint, representing the paint coatings of the ship 
surface. These had also been exposed under the same environmental conditions for the 
same time period as the whole ship surface. 
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TABLE 6 



INSERTED METEOROLOGICAL DATA AND ATMOSPHERIC 

TR.ANSMITTANCE 



DATE 


WIND SPEED 


WIND DIR. 


ICSTL 


ATMOSP. TR-ANSM 


MAY 14 (09:00-09:25) 


4.80 m s 


290 - 300 


**1 

J 


0.8677 


MAY IS (11:57-13:07) 


6.26 m, s 


270 - 280 


3 


0.8539 


MAY 21 (00:05-00:13) 


5.20 m s 


270 - 300 


j 


0.8640 


.MAY 21 (17:18-17:37) 


10.20 m s 


320 


2 


0.8178 


MAY 22 (00:48-01:03) 


9.25 m s 


290 - 300 


j 


0.8282 



Assimiing that the eniissivity was constant over the spectral region involved, 
we used Equation 5.5 derived in Chapter five, for our calculations 



* 0 ■ ■ ^a^^atm 



^0 = 



V^o - 



Making some more assumptions for the laboratoiy environment in which we were 
working, this equation was simplified further. First we had to assume that the ambient 
temperature and the atmospheric temperature were the same in our case. This meant 
that the equivalent source thermal value for the surroundings ( ly ) had to be the 
same as the atmospheric path ( )• So the above equation was then reduced to 



^ 0 ■ ^atm 



Vo - >a) 



Second, for the distance of about 2 m between the detector and specimens being 
measured, we assumed no atmospheric attenuation, so Ty = 1, and the emissivity 
equation was reduced again to 



o 

o 



0 




D - I 



a 



68 



where I'q is the Thermovision 7S0 measured thermal value, is the calibrated 
thermal value for the known ambient temperature and is the calibrated thermal 

value for the known object temperature. 

With this last equation we could calculate the emissivity of the specimens, 
accurately enough, if we knew exactly the ambient temperature and the temperature of 
the specimens. The procedure was started by heating the specimens as unifonnly as we 
could, using a blower, to a temperature giving good thermal images above the ambient 
temperature. This temperature was checked using a thermocouple. With the best 
thermal image on the screen, the temperature of the specimen, the ambient temperature 
and the Thermovision thermal value measurement were recorded. The ambient 
temperature was measured using a conunon thermometer. 

Knowing the e.xact temperatures, the calibration curves with an appropriate 
correction, discussed in Chapter six , were used to get the thermal values ( ) for the 

specimen and ( ) for the surroundings. Therefore, we had what we wanted for our 

calculation. 

The above procedure was repeated several times at different temperatures for 
each specimen. The results were averaged to the following values. 

a. 0.943 first specimen 

b. 0.958 second specimen 

These values are very close to the values we found in different books for the 
emissivity of all kinds of paint, which was 0.95. For the calculations done in this work 
we considered the emissivity of the ship surface to be 0.95, taking the av'erage of the 
two values we found. 

3. Atmospheric Air Temperature 

The atmospheric air temperature measurements were taken from logged 
meteorological data recorded on the ship. The air mass temperatures measured close to 
our station had strong continental influence, so that they were higher than the 
temperatures measured at the ship location and not representative for our situation. 

The atmospheric air temperatures are shown for each date in the data 
presentation tables. 

4. Ambient Temperature 

These temperatures were measured using a conunon mercury' thermometer 
near the open window, where our scanner was operating. Their value was between the 
higher room temperature and the lower outside temperature. 
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The ambient temperature for each date is shown in the data presentation 

tables. 

C. TEMPERATURE DATA EVALUATION AND PRESENTATION 

Inserting the correct parameters, the computer subprogram Evaluaie image was 
used to evaluate the temperature distribution of the ship, for each thermal picture 
recalled from storage. This temperature distribution was presented on the computer 
screen with the standard 8 color scaling combination. 

For better presentation of these temperatures we divided the ship image for each 
aspect in small squares and we tried to assign an average temperature to each square. 
A cartesian system was used to present these temperatures for the starboard side, the 
port side and the stern. These temperatures were averaged for each aspect, for each 
date, in order to have the final temperature distribution for the particular aspect at the 
time period our data were averaged. 

The problem was that the aspect shown by the ship was not always the same, 
because she was drifting and had to keep in station using the engines. So pictures, say 
of the starboard side, taken under several time periods had different aspect angles with 
the line of sight. This caused reflection problems, with thermal radiation from hot 
surfaces such as the stack reflecting from cooler surfaces. Therefore, our scanner 
measured not only the einitted thermal radiation but both the refected and the emitted 
thermal radiation from that surface. So the temperature measurement we were 
evaluating was higher than the correct temperature corresponding to that surface. To 
minimize this problem we tried to average the temperatures only from those pictures 
showing the most perpendicular aspect to the line of sight. 

However the temperatures obtained were not correct, as discussed in detail in 
Chapter six, since the computer-calculated temperatures had still to be corrected using 
the empirical relation already derived. Hence the average temperature for each square 
had to be corrected using the best fit straight line equation corresponding to the 
thermal range in which the data were collected. 

The corrected average temperature distribution for the starboard side, the port 
side and the stern, available for each date, are given in the following Tables 7 to 11. 
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TABLIi 10 
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TABLE 11 

R/V "I’OIN T SUR" TEMPERA rU RE MliASUREMEN l S 
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D. TEMPERATURE DATA ACCURACY 



As mentioned before, six thermocouples were installed about the superstucture, 
at the following locations; 

T1 ; On the after starboard side of the deckhouse 

T2 : On the forward starboard side of the deckhouse 

T3 : On the aft bulkhead of the bridge { away from the stack ) 

T4 : On the aft bulkhead of the deckhouse 

T5 : On the port side of the deckhouse 

T6 ; On the after deckhouse forward of the stack 

These thermocouples were sampled at 10-seconds intervals and the output 
voltages were recorded and later converted into temperatures. Thus, we had an 
indication of the actual temperature at six dilTerent locations on the ship, for the whole 
time period we were taking our measurements. In this way we had an absolute 
reference to compare our data with to estimate the accuracy of the process of 
calculating the temperature of the target from measurements of thermal radiation. 

First we had to fit the thermocouples in the squares covering the areas where 
they were located. Then we had to find the mean calculated temperature corresponding 
to these locations. For e.xample for a thermocouple located at the intersection of two 
or four squares, the calculated temperature corresponding to this thermocouple was 
taken as the mean temperature of the two or four squares. Finally we had to compare 
the above temperature with the thermocouple temperature averaged over for the same 
period as the calculated temperature for each square. 

In the following Table 12 are shown the thermocouple locations and the 
comparison of the actual with the measured and calculated temperatures in degrees C 
for each date. 

As we can see, for most of the time the temperature dilFerences were small. The 
big negative dilTerences in the stern area for the two first dates ( temperature of 
thermocouple higher than AGA thermovision temperature ). are believed to occur 
because the thermocouples were not yet stabilized at that period. That is, at a time 
when the ship was turning, hot air from the stack area was blown onto these 
thermocouples, making their measurements fiuetuate. So for measurements at this time, 
just after the turn, the thermocouples did not have the time to be stabilized. 
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STERN PORT STARE STARE STERN PORT STARE 



table 12 

THERMOCOUPLE - AGA THERMOVISION TEMPERATURE COMPARISON 
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14. 60 


15. 15 


- 0. 55 


9. 53 


10. 10 


* 0. 57 



E. RADIANCE CALCULATION 



Knowing the temperature distribution of the ship, we could calculate the radiance 
for the spectral region of the measurements. 

Starting with the Planck blackbody radiation law we could write 

X, X, . 

/ / 2nhc"^ 1 W 

M = J W(X,T)dX = J ( ^ ) ( ) dX • 

X, X. cm- 



If we let Cj 
follows 



2k1ic- 



and C-) 



he 

we could write the above equation as 

K 



^1 ^1 W 

M = f : ; dX ; 

X^ (e^2 - 1) cm' 



Solving this integral in the spectral region 8-14 ^m. for the known temperatures 
we found the in-band radiant llux density M in W, cm^. With the assumption that our 
surface was an opaque Lambertian surface, we could divide the radiant flux by n and 
multiply by the emissivity ( c ) to obtain the required spectral radiance. 

cM W 

L = . 

■> 

Jt cm"sr 

The integration was done numerically using a program for the T1 - 59 calculator, 
details of which are given in Appendix B. The results of our calculations are shown in 
the following Tables 13 to 17. 
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TAHITI 14 

R/V "POINT SUR" RADIANCIi iVlBASlJRliMIINTS 
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TABLE 15 
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TABLE 16 

K./V "POINT SUR." ITADIANCE MEASUREMENTS 
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VIII. DATA ANALYSIS 



A, INTRODUCTION 

In this chapter we will analyze how the sea background radiance was affected due 
to reflected sky thermal radiation. We will discuss how the reflectance and the 
emissivity of a sea surface roughened by various wind speeds affected the background 
radiance. Also we will calculate the apparent temperature of the sea surface, and will 
find the effective target to background temperature difference sensed by the 
Thermovision under the specific environmental conditions e.xperienced. 

B, EFFECTIVE SEA SURFACE TEMPERATURE 

1. Sea Surface Radiance 

As mentioned in Chapter III the infrared radiance of the sea is the sum of the 
radiance due to thermal emission and the radiance due to reflected incident radiation. 

The reflected incident radiation in the spectral region of 8 - 14 pm is only 
thermal sky radiation by day and night. The contribution of this radiation to the total 
sea radiance depends on the reflectance of the sea surface (R). 

Because the water is essentially opaque in the above spectral region the 
reflectance also controls the emissivity of the sea surface, since for an opaque material 
we can write e = 1 - R . 

Therefore, we can write that the total radiance of the sea surfitce is 

Lj = ( 1 - R ) + R , (eqn 8.1) 

where is the radiance of the sea surface. is the radiance due to sea water thermal 
emission and is the radiance of the sky. To write this equation we assumed that 
the emissivity of the sky, as we have previously discussed in detail, is one. 

2. Rellectance of the Sea Surface 

An important factor in determining the reflectance of a surface is its geometry. 
For the sea surface this geometry depends on the slope distribution of the waves at 
various wind speeds. 

In our discussion about the geometry of the sea surface we have shown, in Fig 
3.8, the difference in reflectance between a flat sea surface and a surface roughened by 
a Beaufort 4 wind speed. 
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These results, as Fig 3.8, were obtained by Cox and Mtink, working in the 
visible range with solar radiation. Such information is lacking for the infrared and night 
conditions, although the above curves seem to be helpful in this area also. We make 
this assumption since the reflectance of the sea water at zero angle of incidence is 
almost equal to the average rellectance in the infrared region from 2-15 pm ( see Fig 

3.4 ), and the sky radiance variation is almost the same for day and night conditions. 

For our analysis we used the same curves corrected for the average index of 
refraction in the spectral region S - 14 pm. Actually we wanted to find a correction 
factor for the reflectance. Since we knew the index of refraction of water for both the 
visible and 8 - 14 pm regions, we could find the corresponding reflectances at normal 
incidence using the equation R(0) = { n - 1 )"'{ n + 1 )^. Their ratio gave us the 
required correction factor. This factor was used to correct the given curves, assuming 
that the correction we found for normal incidence applied for all incidence angles. 

[Ref 7] gives the index of refraction used to calculate the rellectance. for 
uncontaminated sea water, as n = 1.338 and quotes the corresponding calculated 
value of Ry(0) as 0.020. Even though the reflectance equation for n = 1.338 gives 
R.^.(0) = 0.02089, we will use for our calculations the valueof 0.020 used to draw the 
curves. The average index of refraction for the region 8 - 14 pm was found from Fig 

3.5 to be n = 1.303 ; with this index of refraction the reflectance was calculated as 
R-j.{0) = 0.01731. As we can see the reflectance for the region 8 - 14 pm at normal 
incidence is lower than the reflectance in the visible region, in the ratio of Rjj.(0) R.^.(0) 
= 0.865. With this correction the curves were shifted a little lower. 

The following Tables 19 to 23 at the end of this chapter show the reflectance 
(R) values and the emissivity ( 1 - R ) values as they were found from the corrected 
curves, for a flat sea surface and one roughened by a Beaufort 4 wind, for certain 
angles of incidence. Here we reiterate that the values obtained for incidence angles near 
the horizon are less accurate, because the curves are not strictly applicable in this area, 
as discussed in Chapter 111. 

3. Apparent Temperature Calculation 

As we said before in section one, the sea surface radiance can be written as 
the following sum ( Eqn 8. 1 ) 

Ls - ( 1 - R ) + R Lj].y , 



85 



that is, that the radiant flux of the sea surface is the sum of the emitted radiant flux 
and the reflected incident radiant flux. For the whole spectral region we can write 

= ( 1 - R )(7T,/ + R^Tsky'^ ’ 

where is the apparent sea surface temperature, is the actual sea surface 
temperature and is the temperature of the atmospheric air. Dividing both sides of 
the equation by <r we get the more simplified form 

= ( 1 - R ) T,/ 4 - R (eqn 8.3) . 

With this equation we could calculate the apparent sea surface temperature, 
the temperature that our thermal system would sense if it were operating in the whole 
spectrum region and the atmospheric attenuation was neglected. 

With the above assumption, we calculated with moderate precision, for the 
region of observation, the apparent sea surface temperature for a flat sea surface and 
one roughened by Beaufort 4 wind. Calculations were done for incident angles at 89°. 
where we had been working. 85°. 80°, 75° and 70°. The reflectance (R) for the two 
kinds of surfaces and the incidence angles listed was found using the curves in Fig 3.8. 
after being corrected for the spectral region 8 - 14f.im. The actual sea and air 
temperatures were taken from the recorded meteorological data on the ship. 

The results of these calculations are shown in the Tables 19 to 23 at the end 
of the chapter. Detailed analysis and discussion of these results is given in the following 
section. 

4. Calculated Temperature Data Discussion and Presentation 

From the equation we used to calculate the apparent sea surface temperature 
we can see that the two terms at the right side are inversely related to the reflectance of 
the sea surface. 

Therefore, when the reflectance of the sea surface is high, the term with the 
sky temperature ( T^j.y ) dominates, and conversely when the reflectance is low the 
term with the actual sea surface temperature { T^^, ) dominates. So whenever the 
reflectance was high the apparent sea surface temperature we had calculated was closer 
to the sky temperature. Since for our case the sky temperature was always lower than 
the actual sea surface temperature, the apparent sea surface temperatures we have 
calculated were lower than the actual. On the contrarv', whenever the reflectance was 
low, which means that the emissivity of the sea surface was high, the apparent sea 
surface temperatures we have calculated were close to the actual sea surface 
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temperatures, but still some tenths of a degree lower, because of the small contribution 
from the second term. 

Looking at the reflectance values given for the two kinds of surfaces, at 
various incident angles, in Tables 19 to 23 we see that, for large incidence angles (i.e.. 
grading incidence) the difference in reflectance of the two surfaces is large. But as the 
incidence angle decreases the difference decreases and at about 70° goes to zero, {for 
more details see Fig 3.S). For this reason the difference between the apparent 
temperatures we have calculated is large at high incidence angles ( 89°, 85° ) but goes 
to a few tenths at 70 °. 

The calculated temperatures are shown in the T^ row of the Tables 19 to 23. 
In the T^ row' the diflerence between the actual sea surface w'ater temperature and the 
calculated sea surface temperature is shown. T^^.^ is the temperature computed by the 
Thermovision 780 system, divided by the atmospheric attenuation. This was done 
because we wanted to compare this temperature with the T^ temperatures. But since 
the Tj temperatures do not include any atmospheric attenuation, it was necessary 
either to include the atmospheric attenuation into these values or amplify the 
Thermovision measurement, to obtain a value without the atmospheric attenuation. 
We followed the second method and we got the values T^^ shown in the Tables 19 to 
23. 

Comparing these T^^ values with the T^ values we have calculated for incident 
angle at 89°, we can see clearly how much the reflectance of the sea surface, 
roughened at various wind speeds, was affecting the temperature sensed by the 
Thermovision. 

Finally in the row' of the Tables 19 to 23. the difference between the 
actual sea surface temperature and the Tj^.^ temperature is shown. .As we have expected 
this difference is big at low' wind speeds when the reflectance is high and the apparent 
temperature is closer to the sky temperature. 

Generally we can say- that the apparent sea surface temperatures found were 
lower than the actual sea surface temperatures, since the sky was always colder. Their 
difference was related to the incidence angle and the roughness of the sea surface, with 
higher values at incidence angles near the horizon and lower wind speeds. 
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C. TARGET-TO-BACKGROUND TEMPERATURE DIFFERENCE 



1. Effective Temperature Difference 

An important factor for the target detection probability of an infrared system, 
is the target-to-background temperature dilTerence. The higher the temperature 
difference between a target and its background the greater the probability and the 
distance of detection. 

But as we said before the effective temperature difference ( AT ) that our 
thermal system would sense, under various environmental conditions, looking at 
various incidence angles, might not be the same as the actual temperature difference. 

Therefore, the precalculated detection distance, using the actual temperatures, 
at which we expect our system to detect a specific target, may not coinside with the 
actual distance at which the target is detected by the system. 

In this analysis we attempted to find the effective AT, sensed by the 
Thermovision, under the various environmental conditions worked in. and compare it 
with the actual AT . The effective AT was found by taking the difference between the 
temperatures the system calculated for the hottest spot on the ship, ie. the stack, and 
the sea surface. These temperatures, were calculated under the ideal conditions in 
which a thermal imaging system operates, considering emissivity one ( e = 1 ) and 
neglecting atmospheric attenuation ( = I ), and were corrected using the calibration 

relation we had derived. The actual AT was found using the temperature of the stack, 
computed from the radiation measurements using the exact parameters (emissivity. 
atmospheric attenuation, etc. ). 

These temperatures, after being corrected for calibration purposes, were the 
most accurate available for our case, since we did not have any thermocouples on the 
stack to give a direct indication of the actual temperature. The temperature of the sea 
surface was taken from the meteorological data recorded on the ship. 

From this analysis we found that the effective AT, for the environmental 
conditions of the ^measurements and incidence angles near the horizon, was always 
higher, by about a factor of 1.6 to 2 .than the actual AT. This was an indication that 
variations in emissivity and reflectance of the sea surface under the experienced 
environmental conditions had affected the effective AT. sensed by the Thermovision. 

Our results are shown in Tables 19 to 23 at the end of the chapter. A more 
detailed discussion of these results is given in the following section. 
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2. Temperature Difference Data Discussion and Presentation 

The AGA temperatures shown in the Tables 19 to 23, for the hottest spot on 
the ship and the sea surface are temperatures calculated by the computer program 
using the thermal values measured with the Thermovision 7S0, assuming £ = 1 and t 
= 1 and corrected for calibration. 

Therefore, these temperatures correspond to the thermal radiation attenuated 
by the atmospheric elfects and the emissivity of the measured surface. They are the 
effective temperatures our system was sensing under the environmental conditions 
experienced on each date, at high incidence angles near the horizon ( about S9° ). 

The values we present in the Tables 19 to 23 for both the effective and actual 
temperatures are the average values from the data we obtained for each observation 
date. 

Comparing the effective with the actual AT values we got for each observation 
period, we see that the effective AT for the environmental conditions experienced was 
always higher than the actual. It is obvious that there are two main factors causing the 
higher effective AT, the entissivity difference between the two surfaces and the reflected 
incident radiation. Since the emissivity difference was small, the dominant factor was 
the reflected incident radiation. Because the reflectance of the ship surface was too 
small (C = 0.95) to cause such a difference, it was the reflectance of the sea surface we 
had to check. But we have already done this in the previous sections the calculation of 
the apparent sea surface temperature. Then the conclusion was reached that the 
apparent sea surface temperature was lower than the actual sea surface temperature 
because the sky was always colder than the sea for the local conditions. These lower 
temperatures were related to the incidence angle and the roughness of the sea surface, 
with the lower values noticed at incident angles near the horizon and low wind speeds. 

Hence the Thermovision was sensing lower sea surface temperatures due to 
reflected thermal radiation from the colder sky. This was causing the effective AT to be 
higher than the actual. Therefore, we would expect these lower sea surface 
temperatures and high ATs to be related to the wind speed conditions which control 
the roughness and consequently the reflectance of the sea surface, since the incidence 
angle was effectively constant (near the horizon, about 89° ). 

From the tabulated and plotted (Fig 8.1 at the end of the chapter) ratio of the 
elective and the actual AT for the various wind speed conditions we experienced, a 
significant trend can be seen. Even with the assumptions and cumulative errors 
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inherent in the approach it can be seen that for the lower wind speed of 4.S m s on the 
14^^ of May we had an effective AT higher than the actual by about a factor of 2.01 
and for the higher wind speed 10.2 m s on the 21^'- of May ( 17:18 - 17:37 ) we had an 
elTective AT higher only by about a factor 1.63. 

Thus as expected for the lower wind speed, higher reflectance, we had the 
higher AT and for the higher wind speed, lower reflectance, lower AT. It was obvious 
that the effective AT was strongly dependent on the reflected incident sky radiation. 
Also within the spread of the plotted points we may conclude that the AT ratios are 
proportional to the wind speeds. 

D. CALCULATION EXAMPLE 
1. Introduction 

To help understanding of the following Tables 19 to 23. we give here a 
calculation example to show how the data presented were obtained. In this example we 
use the information we had for the IS’’^^ of .May experiment. 

DATE - TIME ; 18^^^ .M.AY 1 1:57 - 13:07 

AIR TEMPER.ATLRE (T^j,^.) : 9.7°C 

SEA SURFACE TEMPER.ATURE (T^^,) : 13.1 '*C 

WIND SPEED - DIRECTION : 6.26 m s(B-4) 280® 

OVERCAST, HAZY, VISIBILITY 3 NM 

ATMOSPHERIC TR.ANSMITTANCE : 0.8539 

Average temperatures computed from Thermovision values for the hottest 
point on the ship (Hot Spot) and the sea surface temperature (S.S.T.), (assuming c = 
1 and = 1) are as follows: 

UNCALIBRATED VALUES : Hot Spot : 22.80°C 

S.S.T : 20.65®C 

CALIBRATED VALUES : Hot Spot : 11.72°C 

S.S.T : 9.52°C . 

All data for this date were collected using THER.MAL lUANGE 2 so the calibrated 
values were found using the straight line ec| nation Y = (1.0238) X - 1 1.62 . 
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The relation between the Beaufort sea state scale and wind speed (nisj is 
shown below in Table IS. 



TABLE IS 

BEAUFORT SCALE 

B WIND SPEED m s 

0 0.00 - 0.46 

1 0.51 - 1.54 

2 2.06 - 3.0S 

3 3.60 - 5.14 

4 5.65 - S.23 

5 S.74 - lO.SO 

6 11.31- 13. S9 



2. Apparent Sea Surface Temperature 

To calculate the sea surface apparent temperature we used Equation 8.3 

T,-* = ( I - R I V + R . 

Using the values in Fig 3.S for incidence angle S9°, where we have been 
working, we got the reflectance values for the flat sea surface and that roughened by a 
Beaufort 4 wind. These values, corrected for the spectral region S - 14 pm using the 
correction factor R- , R.. = 0.S65, are the followins; 

11 V 

FLAT SURFACE R = 0.S04 
ROUGHENED SURFACE R = 0.110 . 

Therefore, from Equation S.3 we get for the flat surface 

Tj^ = { 1 - 0.804 )( 13.1 + 0.110 ( 9.7 )“* 

Tj = 10.65 °C 

and for the rough sea surface at Beaufort 4 wind speed we get 

T^-^ = ( 1 - 0.110 )( 13.1 )“* + 0.110 ( 9.7 )■* 

Tj = 12.S4°C. 

The T^ values represent the dilTcrence ( T^^. - ) between the actual sea 

surface temperature and the apparent temperature we calculated. Hence, we have for 
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the flat surface 13.1 - 10.65 = 2.45°C and 13.1 - 12. S4 = 0.26°C for the roughened 
surface. 

The value is the sea surface temperature computed by the AGA system, 
but without the atmospheric attenuation. To get this value we proceeded as follows. 
The uncalibrated value which the system calculated for the sea surface, 20.65 °C. is 
amplified by the atmospheric transmission = 0.S539; hence, we have 

20.65 4 0.S539 = 24.40°C . 

This value was corrected for calibration purposes using the straight line equation for 
thermal range 2 . to give 

Tj^ = 24.40 ( 1.023S ) - 1 1.62 = 12.S3°C . 

As vve would expect, this value must be close to the T^ value for the roughened surface 
{ B-4 ), since the wind speed on this date was 6.26 m, s ( B-4 ). 

The Tj^^ values are the differences T^^, - T^^, and indicate how much colder 
the effective temperature sensed by the Thermovision was than the actual sea surface 
temperature. Hence we have 

T^j = 13.10 - 12.S3 = 0.27 °C . 

3. Effective Sea Surface Temperature 

As we have said before, the values given for the .AGA are the values the 
system calculated, corrected for calibration . These values are not corrected for the 
atmospheric attenuation and emissivity. Therefore, for the hot spot we have the 
uncalibrated value 22.S°C corrected to 1 1.72°C, using 

AGA HOT SPOT = 22.8 ( 1.0238)- 1 1.62 = 1 1.72°C . 

Similarly, for the value of the AG.A sea surface temperature we have the uncalibrated 
value 20.65 °C which corrected gives the value 9.52 °C. 

For the actual conditions, the temperature of the hottest spot on the ship, 
since we did not have any thermocouples on this location . was taken as the highest 
temperature we found for this location in the calculated temperature distribution of the 
ship. This value for the 18^^^ of May was 14.28°C. The sea surface temperature was 
taken from the recorded data on the ship. 
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The AT value represents the hot spot to sea surface temperature dilTerence. 
Hence, we have effective AT^ = 2.20 °C. from 

ATg = 1 1.72 - 9.52 = 2.20°C . 

The actual AT^^ = 1.18°C, from 

AT^ = 14.28 - 13.10 = 1.18°C . 

The ratio ATg AT^^ = 1.864 represents the fractional increase in elTective temperature 
for the particular environmental conditions encounered. 
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TABLE 19 

SEA SURFACE TEMPEFL-\TURE AND EFFECTIVE AT "A" 



DATS / TIMS 

WIND SPSSD / DIRSCTION 
AIK TSMKSKATURS 



14 MAY 09: 00 - 09: 25 
4.8 m/sec (3-3) 290»- 300» 
10. S»C 



AMBISNT TS.MPSRATUKS ; 15 »C 

SSA 3'JKSACS TSMKSKATJKE : 12. 6»C 
OVSKCA.ST, HAZY, VISIBILITY 4 NM 
ATMOSKKSRIC TKANSMITT.ANCS : 0. 3679 





ANGLE 


OF INCIDENC: 


: AND SURFACE CONDITIONS 


• i 


i 8S° 




1 so" 


j 75° 1 70° i 


' 1 


(FLAT : 5-4 


1 I FLAT B-4 


' 1 FLAT ' 5-4 


MfLAT 1 3-4 ' FLAT 3-4 


' R 1 


lo. 80410. 110 


i lo. 513 IC. 1381 


i 0. 292 i 0. 166 


'10.18010.151^0.120:0.120 i 


1-R j 


0. 196|0. 390 


1 0. 487 1 0. 862 , 


iC. 70810. 834 


jo. 32010. 3491 !0. 83010. 330 ' 


1 1 

Ts i 


11. 01 1 12. 42 


ill. 66ll2.36^ 


’ 12. 09 j 12. 32 ! 


' 12.29112. 35I il2. 40! 12. 40 ; 


: Td 1 


' 1. 59 ' C. 18 


i 0.94 0.24' 


0.51 C.2S^ 


0.31| 0. 2sj i 0. 20| 0. 20 ' 



Tm 


1 11. 51 i 1 


T^d 


1 1. 09 ! 



1 TEMPER. 1 


HOT SPOT 


1 SSA SURFACS 


AT 


i RATIO ' 


' AGA ! 


IS. S2 


1 3. 44 


7. 38 


I 2. 01 


ACTUAL i 


16. 27 


! 12.60 i 


3. 67 





R : RSSLSCTANCS OS TKS SSA SURFACS 

FLAT: FLAT SSA SURFACS 

3-4 : SSA SURFACS ROUGKSNSD 3Y BSAUFORT 4 WIND 

Ts : GALCULATSD APPARSNT SSA SURFACS TSMPSRATURS 

Td : DIFFSRSNCS 3STWSSN ACTUAL AND CALCULATSD SSA TS.MPSRATURSS 

Tni : AGA StA TSMPSRATURS MSASUREMSNT WITHOUT ATMOSPHERIC ATTENUATION 

Tmd : DIFFSRSNCS BETWEEN ACTUAL AND AGA SEA TEMPERATURES 



* ALL TEMPERATURES ARE IN °C 
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TABLE 20 

SEA SURFACE TEMPER.ATURE AND EFFECTIVE AT "B" 



/ TIME 



: 13 11:57 - 13:07 



WIND SrEZD / DIEECDICN 



5.25 n:/sec (3-4) 270o-280o 



AIR TE:*!?ERArjRE 



9. 7»C 



;^2^3:EMT TEr-IrERArjRE : 17oC 

SEA SUREACE TEriRER-i.rjRE: : 13.1»C 

0^v.^RCAST, HA2Y, VISIBILI’TY 3 NW 



.4D:^.C SPHERIC TRANSHirC.^^TCE : 0.3539 



1 

1 


.ANGLE 


0? IMCITHTJCT SURFACE CONDITIONS 


1 


‘ [ 39’ 


' ! 35’ 30’ 1 i 75’ ! i 70’ 


, - 


’ FLAT 1 3-4 


: ELAE 3-4 ELAE ' 3-4 ’ lELAE ! 5-4 ’ ! ELAE 3-4 


R 


'' 0. 304! 0, 110 i 


'0. 31310. 133' 0. 292 !0. 156m 0. 13010. 131: lo. 120 iC. 120 



1-R ' 0. 196 0. 390 ■ ' C. 487 ' 0. 352 » ' C. 708 ' 0. 334 i ! 0. 320 10. 349 i ! 0. 38010, 330 



Es 


* 10. 35 1 12. 341 


' 11. 72|12. 771 


i 12. 37 j 12. 70 i 


i 12. 561 12. 74j 


1 12. 32 1 12. 32 


Ed 


1 2.45. 0.25| 


: 1.38i 0.331 


! 0.73i 0.40 1 


1 0.441 0.3oi 


!• 0. 2Si 0. 23 i 



i Er\ 


12. 33 1 i 


|Tnd 


0.27 h 



EEMEERv. 


HOE SPOE 


SEA SURE.ACE 


AE 


RAEIO 


.AG. A 


' " . 72 


9. 52 


2. 20 


1. 36 


ACEUAL 


14. 23 


12. 10 1 


1. 13 





R : REEIECD.ANCE OE THE SEA 3URE.-.CE 
ELAT: EDAE SEA SUREACE 

3-4 : SEA SURE.nCE ROUGHENED 3Y 3E.VJEORT 4 WIND 

Es : CAECYIAEED ARP.^^RENE SEA SUREACE EE:^ER.AEU?vE 

Ed : DIEEERENCE SETWEEN ACEUAL .A2ID GALCULAEED SEA EENPER.AEURE3 

E.*n : .AG A SEA EENPERAEURE MEASUREMENE WIEHOUE .AENO SPHERIC .AEEElfJAEION 

Er.d : DIEEERENCE 3EEWEEN .ACEUAD .AND AGA SEA EEMPERAEURES 



.ALE EENPERA.EURES ARE IN ®C 
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TABLE 21 

SEA SURFACE TEMPER.ATURE AND EFFECTIVE AT "C" 



DATE /TIME : 

WIND SPEED / DIRECTION : 
AIR TEMPERATURE : 

AMB I ENT TE:^ ERA.TURE 
SEA SURFACE TEMPERATURE : 
CLEAR, VISIBILITY 9 IJM 
ATMOSPHERIC TRANSMITTANCE 



21 MAY 00: 05 - 00: 13 

5.2 m/sec (3-4) 270«- 30Qo 

11. 3«C 

17«C 

11. 9°C 

: 0. 3640 



1 


ai:gle 


OF INCIDENCE 


.AND SURF.ACE CONDITIONS 


1 


1 


!i 39’ 


35’ 


90® 


75® ! 


70’ ! 


! i FLAT 


1 3-4 




FLAT ■ 3-4 i 


IFLAT ' 3-4 


' 1 FLAT 


3-4 1 


IfLAT ! 3-4 ' 


1 

R 


' jo. 804i0. 110' 


0. 513 i 0. 138 1 


|0. 292(0. 166 


: 0, ISO 


10. 151j 


0 . 120 jo. 120 


1-R 1 


1 io. 196 


Io. 890 i 1 


1 ' 

0. 487 1 0. 862 1 1 


0. 708i 0. 834! 


1 |o. 820 io. 849 i 1 


0.880(0.880 


Ts ^ 


1 11. 42 


11. 84i 


1 


11. 60 i 11. 82 ! 


11. 73 ! 11. 81 


11. 80| 


11. 81 1 


11. 83(11. 83 j 


Td 1 


1 1 

1 0.48' 


0. 06 1 ' 


0. 30 1 0. 08 1 i 


0. 17 1 0. 09 1 


1 0 . lol 


0. 09 M 


0.07! 0.07 5 



Tm 


1 11. 65 j 


Tmd 




0. 25 





1 i 

TEMPER, i 


HOT SPOT 


SEA SURF.4CE 


AT 


I RATIO 


1 AGA 1 


12. 88 


7. 74 


! 5. 14 


1. 90 ! 


ACTUAL i 


14. 60 


1 11. 90 1 


2. 70 


! 1 



R : REFLECTANCE OF THE SEA SURFACE 

FLAT: FLAT SEA SURFACE 

3-4 : SEA SURFACE ROUGHENED BY BEAUFORT 4 WIND 

Ts : CALCULATED APPARENT SEA SURFACE TEMPERATURE 

Td : DIFFERENCE BETWEEN ACTUAL AND CALCULATED SEA TEMPERATURES 

Tm : AGA SEA TEMPERATURE MEASUREMENT WITHOUT ATMOSPHERIC ATTEIfUATION 

Tmd : DIFFERENCE BETWEEN ACTUAL AND AGA SEA TEMPERATURES 



* ALL TEMPERATURES ARE IN '’C 
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•TABLE 22 

SEA SURFACE TEMPEFL^TURE AND EFFECTIVE AT "D 



DATE / ZVAZ 

WI>ro SPEED / DIRECTION 



21 MAY 17: 18 - 17: 37 
10.2 m/sec ( 3-5) 320<> 



AIR TEMPERATTJRE : 12.2oC 

AMBIENT TEMPERATURE : 19 oc 

SEA SURE.ACE TEMPERATURE : 13.2oC 

SUN, CLEAR, VISIBILITY 10 NM. 



ATMOSPHERIC 



TRANSMITTANCE 



0. 8178 



ANGLE OF INCIDENCE AND SURFACE CONDITIONS 





' * 39° ' ' 35° 30° 




1 75° 


M 70® 


; -IrlAT . 3-4 liiLAT 3-4 j | 3 LAT 1 3-4 


1 


ItLAT j 5-4 


i i FLAT 1 3-4 j 


R 


! 0. 304 i C. :i0| |0. 513 1 0. 133 i i 0. 292 i 0. 166 




0. 180|0. 151 


jo. 12010 . 120 


1-R ' 


io. 196-0. 390 1 lo. 48710. 362 | | 0. 708 1 0. 83 41 


!o. 82010. 849 


lo. 88010. 880 1 


t Ts i 


12. 41 , 13. 10 M 12. 72 1 13. 0711 12. 93 | 13. 05 i 




! 13. 04' 13. 06 


1 13. 09 ! 13. 09 1 


Td j 


! 0.79 1 0.10'i 0.4sl 0. 13 h 0.271 0.15 


! 0. 16' 0. 14| 


! 0. 11 1 0. 11 ' 



1 Tm 1 ’ 


13.13 -I 


1 T.Tid 1 ’ 


0. 02 I i 



TEr^ER. ' 


HOT SPOT 


' SEA SURFACE i 


dT 


i RATIO 1 


AGA 


15. 30 


1 

j 9.26 


5. 04 


! 1. S3 


! ACTUAL 


15. 89 


13. 20 


3. 69 


■ i 



R^ : REFLECTrPICE OF THE SEA SURFACE 

FLAT: FLAT SEA SURE.ACE 

3-4 : SEA SURFACE ROUGHENED 3Y 3EAUF0RT 4 WIITO 

Ts : CALCULATED APPARENT SEA SURE.ACE TEMPERATURE 

Td : DIFFERENCE BETWEEN .ACTUAL .AID CALCULATED SEA TEMPERATURES 

Tni ; AGA SEA TEMPERATURE MEASUREMENT WITHOUT ATMOSPHERIC ATTENUATION 

Tnid : DIFFERENCE 3ETV7EEN ACTUAL AND AGA SEA TEMPERATURES 



.ALL TEMPERATURES ARE IN ®C 
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TABLE 23 

SEA SURFACE TEMPER.ATURE AND EFFECTIVE AT "E" 



DATE / 

wi:td speed / dipection 
AIR te:c?erature 

AT'IBIENT TEr^ERATURE 



22 yiKY 00: 48 - 01: 03 
9. 5 m/sec (3-5) 290° 
9. 7°C 
17°C 





SEA SURFACE 


lEMPERATURE : 


10. 3°C 






CLEAR, VIS 13 


IT ITU 10 NM 








.ATMOSPHERIC ' 


CRANSMITT.ANCE 


: 0. 3232 




i 


ANGLE OF INCIDENCE 


AND SURF.ACF CONDITIONS 1 


■ ' 1 39“ 1 


1 85“ i 


1 80“ i 


1 0 M 0 1 

1 /O i 1 /O 1 


i 


i FLAT 1 3-4 


IFLAT ’ 3-4 1 


(FLAT ! 3-4 1' 


FLAT i 3-4 1 i FLAT ! 3-4 j 


R 1 


io. 30410. lioi 


lO. 313|o. 1331 


io. 292i0. 165i 1 


0. ISOlO. 13ll io. 12010. 120 ' 


1-R 1 


0. 196 i 0. 890 i 


10. 48710. 362 i 


jo. 708 io. 834j 1 


0. 32010. 349 1 !0. 38010. 380 ' 


Ts 1 1 


1 9. 94110. 59 1 


10. 281 10. 57 1 


1 10. 3l! 10. 54i i 


10. 52 i 10. 55 1 1 10. 63 i 10. 53 1 


'^d i i 


0. 36i 0. ui ! 


0.32 0.13 1 


i 0. 29 1 0. loM 


0. 18l 0. 15 1 ! 0. 12 i 0. 12 1 


1 1 

T:n 1 1 


10. 76 i 1 








Tmd 1 


0. 04 1 1 









1 TEi^ER. 1 


HOT SPOT 


! SEA SURFACE j 


AT 


RATIO 


j AGA i 


3. 54 


1 5. 30 1 


1. 74 


1 1. 70 


' .ACTUAL i 


11. 32 


' 10. 30 


1. 02 


i 



R ; REELECT.^iNCE OF THE SEA SURFACE 
FEAT: FEAT SEA SURFACE 

3-4 : SEA Su’RFA 4CE ROUGHEMED 3Y 3EAUF0RT 4 WHfD 

Ts : CALCULATED APPA.RENT SEA SURFACE TEMPERATURE 

Td : DIFFERENCE BETWEEN ACTUAL .AND CALCULATED SEA TEMPERATURES 

T.-n : AG.A'* 5EA TEMPEIEATURE .MEASUREMENT WITHOUT .ATMOSPHERIC .ATTENUATION 

T.T.d : DIFFEP^NCE 3ETV/EEN .ACTUAL AND .AGA SEA TEMPERATURES 



* kLL TirtPcSATURiS ARE IN °C 
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WIND Sl'l';i':i) TO 1)T KATIO KI'IDATION 





' , 1 1 

1 

C ; i; ! , 




^ / 

- !/c 




i 

! 1 

i i 

t 


1 

1 


i 

1 




1 1 / : 
1 


1 


i 


1 

! 

i 


1 

1 




j 








1 

i 


1 
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! 


1 

1 


1 

1 


1 


1 

1 ! 






1 

1 

1 


1 


I 


1 

i 

1 


1 ! j 




1 




; 



in ; 

Oi 

lii 

S>l 

^6' 

iSxi 
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in 



in 
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Figure 8.1 A 1 ratios ploled vcisus wind speed (n\/s). 



IX. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

Temperature and radiance distributions of an instrumented target measured 
remotely by radiometric imaging with the AGA Thermovision 780 system have been 
compared with locally measured environmental parameters. The ship-to-background 
temperature dilTerence, frequently used as a measure of the contrast, was compared as 
measured thermometrically with local temperature sensors and radiometrically using 
the Thermovision. Significant differences have been observed, depending on the time of 
the day, windspeed, and atmospheric conditions. These difterences have been assigned 
to effects of the low emissivity and high reflectivity of the sea surface at near grazing 
angles. 

Reflected thermal radiation from the sky contributed greatly to the sea surface 
radiance sensed by the system. Comparison of target-to-background temperature 
differences from radiance and from thermometric measurements showed large 
discrepancies due to the effect of the reflected sky radiation, characteristic of the 
atmospheric temperature. This reflected radiation was strongly dependent on the 
roughness of the sea surface under conditions of varying wind speed. The differences in 
the air and sea surface temperatures caused the effective (radiometric) temperature 
differences between the target and background to deviate both positively and negatively 
from the true temperature difference. With the sky colder than the sea surface, as was 
alwavs found in the first series of measurements, the reflected skv radiation caused the 
sea surface temperature measured by the Thermovision to be lower than the actual. .As 
a result, the effective temperature difference as observed with the Thermovision was 
higher than the actual difference. Both the temperature decrease and the increase in the 
effective temperature difference were related inversely to the sea state conditions. 

An e.xpression for the sea surface radiance including the reflected sky radiation 
has been developed and used to model the radiometric effective ship-to-background 
temperature difference. This model can be evaluated for known atmospheric conditions 
using the LOWTRAN 6 computer code and has been found to be consistent with 
observations to date. 
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In analysis of radiometric images of the ship taken at difTerent aspect angles, 
apparent differences in radiance distribution were observed. This is considered to result 
from thermal radiation from hotter portions of the ship reflecting from cooler regions, 
increasing the apparent surface temperature. Consequently, temperature values from 
only those surfaces more nearly perpendicular to the line of sight were averaged to 
obtain representative temperature and radiance distributions. 

Although the quantitative results cannot be considered to be of great accuracy, 
they show clearly the strong dependence of the reflected sky thermal radiation on the 
sea roughness due to the wind speed at the times of the observations. 

B. RECOMMENDATIONS 

Since we had to calibrate our system empirically to achieve the accuracy of our 
measurements it is obvious that we have inserted some error factors into our results. 
Therefore, we recommend a new series of target and background measurements to be 
taken with the whole system being accurately calibrated. 

A few months before this work was completed the Thermovision was tested and 
calibrated accurately. The gain of the preamplifier in the scanner was adjusted and new 
calibration curves were derived. However there remained a problem in that the 
computer did not calculate the exact temperatures because the thermal values inserted 
by the system into the computer program were off by four thermal units, even though 
a few weeks before this work was completed we took a new series of measurements of 
the same target. Our results and analysis are shown in .Appendix .A. 

.All our measurements were taken at an incidence angle near the horizon where 
the curves in Fig 3.S. used to find the reflectance, are not strictly applicable. Therefore, 
depending on the operational use, we may need this information with greater accuracy, 
we recommend: 

For systems like missiles or airborne forward looking infrared devices a series of 
measurements should be taken for incidence angles between 70 ° and 80° where the 
curves in Fig 3.8 are applicable, and where more accurate results can be obtained. 

For forward looking devices operating on ships, with incidence angles near the 
horizon, long series of measurements, under all possible kinds of environmental 
conditions, may be needed to provide sufficiently accurate information. 
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APPENDIX A 

DATA ANALYSIS OF LATER MEASUREMENTS 



1. GENERAL 

A second set of thermal imaging data were collected in rnid November, 1986. All 
measurements were taken in the spectral region 8 -14 pm as before, using the AGA 
Thermovision supplemented with a microcomputer. The ship target was again the R V 
"POINT SLR", scheduled on an operational oceanography student cruise, and the 
geographic location where the thermal imaging device and the ship were positioned 
during the e.xperiment was also the same. 

These data were collected in order to obtain further, more accurate, 
measurements, under different environmental conditions, after the Thermovision had 
been calibrated. This was done to check our previous results and improve the 
information on the reflected sky thermal radiation and the effective target-to- 
background temperature difference relation. 

2. CALIBRATION TEST 

The Thermovision was accurately calibrated after the first measurements had 
been taken. The preamplifier gain of the scanning device which was found to be 
incorrectly set, was adjusted and new calibration constants were derived and inserted in 
the computer program. However, even if the Thermovision was operating accurately on 
this occasion, the temperatures calculated by the computer were not accurate. The 
problem lay in the thermal values transmitted to the computer from the Thermovision. 
The thermal values received by the computer were about 3 to 4 thermal units lower 
than measured directly on the Thermovision. The error was minimized by adjusting 
the THER.VIAL LEVEL control to read about 4 thermal units higher, using the 
THER.MAL LEVEL ADJ control on the Black, White monitor chassis. In this way the 
the computer received almost exactly the thermal values corresponding to the object 
in view. Of course this was not the right way to solve the problem, since now the 
Thermovision was measuring thermal values 4 thermal units higher than previously. 
However, the system was empirically calibrated again, following exactly the same 
procedure described in Chapter VI. 
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The calibration temperature measurements for the two thermal ranges used in 
these series of measurements { 2 and 5 ) are shown in the following Tables 25 and 26. 
The data were plotted as shown in the Following Figures A.l and A. 2 and the best fit 
straight line was found. .As can be seen from the mean and variance of the dilTerence 
between the e.xact and calculated temperatures, the temperatures the computer 
calculated were ver>' close to the exact temperatures. The correction error found this 
time was within 2.9 " o for the thermal range 2 and 3.5 '!'o for the thermal range 5. Of 
course again we have to keep in mind that these calibration functions are valid only for 
temperatures in the region from 5°C to 25 °C. 

3. MEASURED DATA 

Thermal images of the target and its background were collected on the of 
November ( 09:07 - 09:12 ), on the IG*-^ of November ( 17:58 - 18:02 ) and on the 18’-^^ 
of November (08:34 - 08:36 '). .All measurements were taken using the same equipment 
as before, in the same spectral region and geographic location as is described in 
Chapter Vll. The only difference was that this time the R "POINT SUR"did not 
maneuver near the buoy to show difierent aspects but passed by the buoy on her way 
into and out of the Monterey harbor. The distance between the lab station and the 
ship was about 650 m on the 14^^ and 18’-^^ when the ship was going out of Monterey 
Bay showing the port side and about 750 m on the 16'’^^ on her way into the harbor 
showing the starboard side. 

The environmental conditions experienced during these measurements are shown 
in the data presentation tables for each observation period (Tables 27.28.29). 
Generally we had low speed winds blowing from the North West, overcast skies and 
fog on the last day. The air temperature was always higher on this occasion ( within 
0.5 “C to 1.9‘’C ) than the sea surface temperature. 

To analyze the data the same procedure described in Chapter VI 11 was followed. 
The apparent sea surface temperature was calculated and the efiective target-to- 
background temperature difference, sensed by the Thermovision under the local 
environmental conditions, was found. These results are presented in Tables 27.28,29 at 
the end of this appendix. .Also a plot of the ratios of the efiective to actual ATs versus 
wind speed is presented in Figure A. 3 after the tables. 
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TABLE 24 



THERMAL RANGE 2 

NEW CAL1BR.ATI0N TEMPEFLATURE DATA 



TEMPER.ATURES IN DEGREES C 



BLACKBODY 

EXACT 


AG A 

COMPUTED 


DIFFERENCE 


5.0 


4.7 


0.3 


6.0 


5.8 


0.2 


7.0 


6.7 


0.3 


8.0 


7.8 


0.2 


9.0 


8.7 


0.3 


lO.O 


9.7 


0.3 


11.0 


11.1 


-0.1 


12.0 


12.2 


-0.2 


13.0 


12.8 


0.2 


14.0 


13.7 


0.3 


15.0 


14.9 


0.1 


16.0 


15.9 


0.1 


17.0 


17.2 


-0.2 


18.0 


18.2 


-0.2 


19.0 


18.9 


0.1 


20.0 


20.2 


-0.2 


21.0 


20.9 


0.1 


22.0 


22.0 


0.0 


23.0 


23.0 


0.0 


24.0 


24.1 


-0.1 


25.0 


25.2 


-0.2 



MEAN DIFFERENCE ; 0.02S6 
STANDARD DEVIATION : 0.2305 
VARIANCE : 0.0506 

BEST FIT STRAIGHT LINE ; Y = (0.9798) X - 0.3628 
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[•iguic A. I I hcrm:il Range 2 new calibialion lempeiatine data. 



TABLE 25 



THERMAL R.AXGE 5 

NEW CALIBR.AT10N TEMPER.ATURE DATA 



TEMPER.ATURES IN DEGREES C 



BLACKBODY 

EXACT 


AGA 

COMPETED 


DIFFERENCE 


5.0 


4.6 


0.4 


6.0 


5.5 


0.5 


7.0 


6.6 


0.4 


8.0 


7.6 


0.4 


9.0 


8.4 


0.6 


10.0 


9.6 


0.4 


11.0 


10.5 


0.5 


12.0 


11.4 


0.6 


13.0 


12.5 


0.5 


14.0 


13.6 


0.4 


15.0 


14.6 


0.4 


16.0 


15.7 


0.3 


17.0 


16.7 


0.3 


18.0 


17.7 


0.3 


19.0 


18.6 


0.4 


20.0 


19.6 


0.4 


21.0 


20.7 


0.3 


22.0 


21.7 


0.3 


23.0 


22.6 


0.4 


24.0 


23.7 


0.3 


25.0 


24.6 


0.4 



MEAN DIFFERENCE : 0.4047 
STANDARD DEVIATION : 0.0921 
VARIANCE : O.OOSI 

BEST FIT STR.MGHT LINE : Y = (0.9920) X - 0.5217 
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riguic A. 2 I hermal Range 5 new calibialion leinpeiaUire dat:i 



4. 



DISCUSSION AND ANALYSIS OF THE RESULTS 

As was expected, reflected thermal radiation from the warmer sky affected the sea 
background radiance under the various sea state conditions due to the wind speeds 
encountered. Therefore, the apparent and the efTective sea surface temperatures were 
found to be higher than the actual. As before, the differences between the apparent and 
actual sea surface temperatures were higher at incident angles near the horizontal and 
at lower wind speeds. 

Since the effective sea surface temperature sensed by the Thermovision was 
higher than the actual, the effective AT sensed by the detector was lower than the 
actual AT. 

The major error factor occurred because we did not have an indication of the 
actual temperature of the stack, the hottest spot on the ship, and we had to trust again 
the temperature the system calculated for this spot. On this exercise we did not have 
thermocouples positioned at various locations on the ship as we had before, to 
calibrate the measurements. An attempt was made to record a measure of the stack 
temperature, using a common thermometer fastened to the stack surface. 
Unfortunately since we could not place the thermometer on the surface viewed by the 
system, we did not get accurate measurements. 

Even with all the uncertainties about the exact ATs and their ratios, the results 
show again that the variations in emissivity and reflectance of the sea surface under the 
local environmental conditions had affected the effective AT sensed by the 
Thermovision. 

From the tabulated and plotted ratios (Fig A. 3) of the effective and actual ATs it 
can be seen that for the lower wind speed of 1.5 m s ( Id’’^ Nov.), and the higher 
reflectance, we find the lower AT ratio of 0.939 . This ratio increases to 0.954 when 
the wind speed increases to 2 m. s ( IS'’*^ Nov. ) and the reflectance decreases. Normally 
one would expect to have the higher ratio on the 14^^^ of November when the wind 
speed was 2.9 m s, but unfortunately the measured data do not seem to be consistent 
for that date. The temperature value recorded for the stack on this day appears 
inconsistent with the reliably measured values of air temperature. However, the ratio is 
close to one, not too far from what it is expected. 
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5. CONCLUSIONS AND RECOMMENDATIONS 

It is noticable again in this series of measurements that rehected sky radiation 
afTected the sea surface radiance sensed by our system. Since the sky was warmer than 
the sea surface . during the time periods of observation, reflected sky thermal radiation 
caused the sea surface temperatures sensed by the Thermovision to look higher than 
the actual sea surface temperature. The higher effective sea surface temperatures 
caused the effective ATs sensed by the Thermovision to be lower than the actual ATs 
under the various wind speeds experienced. 

To improve the accuracy of the data set. it is reconunended that a scries of 
measurements be taken, with the whole system correctly calibrated and with a 
measurement of the actual temperature of the side of the stack viewed by the system. 
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TABLE 26 

SEA SURFACE TEMPER.^TURE AND EFFECTIVE AT 'T 



DATE /TIME : 

WIND SPEED / DIRECTION : 
AIR TEMPERATURE : 

ambient: temperature : 

SEA SURFACE TEMPEPA^TURE : 
TARGET DISTANCE 
OVERCAST, VISIBILITY 4 NM 
ATMOSPHERIC TRAIISMITTPN^CE 



14 NOV 09: 07 - 09: 12 

2. 9 m/sec ( 5-2) 320® 

13. 9®C 

15®C 

12.2®C 

550rr. 

: C. 8864 



! 


ANGLE 


OF INCIDENCE AND SURFACE CONDITIONS 


j 


89® j 


! 85® 


!i 8C“ '1 75' 


N 70' 1 


1 1 


IFLAT 


1 B-4 1 


f % m 


i E-4 


! FLAT i B-4 1 1 FLAT ' 5-4 


1 i FLAT : B-4 ! 


I' 1 

1 ^ 1 


0. 804 


0. 110 


1 0. 513 


C. 138 


i 0. 292 10. 166i i 0. 18010. 151 


io. 12010. 120 i 


1-R 


0. 196 


C. 890 


1 C. 487 


i 0. 862 


i 0. 708 1 0. 834 1 0. 820 | 0. 849 ! 


10. 88ojo. 880 


1 Ts 


1 13. 60 1 


12. 34! 


1 13. 11 


12. 29 


I12. 71|12. 45 ! ! 12. 48 ! 12. 42 


i 12. 36i 12. 36 


t Td 1 


; -1. 40 1 


-0. 14 j 


! -0. 91 


' -0. 19 1 


j-0.5l!-0.25 '-O.2sl-O.22j 


1 -0. I61-C. 16 


Tm 


i 


56 


1 








i 1 

! Tmd 1 


1 -0. 


36 











1 TEMPEPs. 


! HOT SPOT 


j SEA SURFACE 


1 AT 


RATIO 


AGA 


15. 00 


11. 13 1 


3.87 


1. 04 


i ACTUAL 


15. 90 


12.20 


1 2.70 1 


1 



R : REFLECTANCE OF TriE SEA SURFACE 
FLAT: FLAT SEA SURFACE 

5-4 ; SEA SURFACE ROUGHENED EY BEAUFORT 4 HINT) 

Ts : CALCULATED APPARENT SEA SURFACE TEMPERATURE 

Td : DIFFERENCE BETWEEN ACTUAL AND CALCULATED SEA TEMPERATURES 

Tm : AGA SEA TEMPERA.TURE MEASUREMENT V7ITHOUT ATMOSPHERIC ATTENUATION 

Tmd : DIFFERENCE BETWEEN ACTUAL AND AGA SEA TEMPERATURES 



* ALL TEMPERATURES ARE IN 'C 
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TABLE 27 

SEA SURFACE TEMPER^MURE AND EFFECTIVE AT "G" 



DArZ /TIME : 

speed / direction •. 

AIR TEMPERATURE 
AMBIENT temperature: : 

SEA SUREACE TEMPERArJPE: : 
T.ARGET DISTANCE : 

PARTIALLY CLOUDY, VI SIBIL 
ATMOSPHERIC TRANSMITTANCE 



15 NOV 17; 58 - IS; 02 

1.5 m/sec (3-1) 32C' 

IS. S'C 

17. 5 ®C 
14'C 
750iti 

ITU' 5 m 
: 0. 3842 



I .ANGLE OF INCIDENCE AND SURFACE CONDITIONS 


ss* * 1 es* i . 80 ® 1 


! 7 


5® 1 , 70® 


j ( ELAT B— 4 


1 1 ELAE : 5-4 


1 ELAE 5-4 i 


1 ELA.T 


! 5-4 


! ELAE , 5-4 


: R 1 


1 C. 8041 C. lie 


i jo. 513 Ic. 1381 10 . 29210. issl 


|0. 180, 


1 0. 151 i 


Ic. 120 IC . 120 1 




ic. 196lo. 890 


| 0 . 487|o . 862 




|o. 70e|o. 834 


0. 820 


0. 849 


io. SCO 10.880 1 


1 Ts 


15. 58114. 25 


1 15. 06 1 14. 31 




14. 63 1 14. 74 


Il4. 40 


14. 39 


i 14. 27 1 14. 27 


' Td 


i -1. S8i -0. 25 1 


1 -1. 06 1 -0. 31 ' ' 


-0. 63 1 -0. 74 I 


! 1 
i -0. 40 i 


- 0 . 39 ! 


l-O. 27 j-O .27 1 



! Tit. | 


i 14. 95 


iTmc 


1 -0.95 



TEMPER. 


EOT SPOT i 


SEA SURFACE i 


AT 


RATIO 


1 


1 17.38 


1 13.22 


4. 16 


0. 939 


i 

1 ACTUAL 1 


1 18.36 


13. 95 1 


4. 43 





F IjAT 



5-4 



Ts 

Tc 



Tmd 



REFLECTANCE OF THE SEA SURFACE 
FLAT SEA SURFACE 

SEA SURFACE ROUGHENED BY BEAUFORT 4 WIND 

CALCULATED APPARENT SEA SURFACE TEMPERA.TURE 

DIFFERENCE BETWEEN ACTUAL AND CALCULATED SEA TE^U’E?RTURES 

AGA SEA TEMPERATURE I'EIASUREMENT WITHOUT ATMOSPHERIC ATTENUATION 

DIFFERENCE BETWEEN ACTUAL AND AGA SEA TEMPER.ATURES 



★ 



TE.'C? ERArj?.£S ARE IN 'C 
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TABLE 28 

SEA SURFACE TEMPER.MURE AND EFFECTIVE AT "H 



DATE / Tir-E 

WIND SPEED / DIRECTION 
AIR TEMPERATURE 
AMBIENT TEMPEFATURE 
SEA SUREACE TEMPERATURE 
TARGET DISTANCE 
OVERCAST, FOGGY, VISIBi: 
ATMOSPHERIC TRANS MITT AN: 



: IS NOV 06: 34 - OS: 35 

: 2.0 m/sec (5-2) 335* 

: 14. S*C 

: 15"C 
: 14. 3S®C 

: 650r. 

ITY 1 NM 
E : 0. 8890 



! 


ANGLE 


OF incident: 


1 SURFACE CONDITIONS 






1 89' 


! 1 es' 


: ' 80® 


75' 


70® 


i 


FLAT 1 


i 5-4 


1 1 FLAT 


5-4 


1 1 FLAT ( 5-4 


' IfLAT i 5-4 


1 1 FIA.T 


5-4 : 


: R 


0. 804 


0. 110 


10. 513 


lo. 138 


1 10.29210. 166 


1 Ic. 180 l 0 . 151 


i 0. 120 


10.120 1 


!i-r 


Ic. 196 


0. 890 


lo. 487 


10. 862 


lo. 708 |o. 834 


lo. 820 lo. 849 


lo. 880 


0. 880 


' Ts ! 


i 14. 72 ! 


14. 42 


I 14. 60 ! 


14. 44 


i 14. sol 14. 4S 


• i 14. 46| 14. 44 


' 14.43 1 


14. 43 i 


i Td 


1 1 

! -0. 34' 


1 

o 

o 


! -0. 22 : 


1 

o 

o 


! -C. 12 ' -C. 07 


1 

o 

0 

CO 

1 

o 

o 

O' 


i -0. 05 i 


-0. 05 j 



Tm 


! 14. 62 


j Tmd 


1 -0.24 ! 



TEMPER. 1 


HOT SPOT 


SEA SURFACE 


AT 


i RA.TI0 


! AGA 


14. 67 


!-• 

o 

o 


1. 67 


1 1 

1 0. 954 1 


ACTUAL 


16. 13 


14. 38 


1. 75 





R : REFLECTANCE OF THE SEA SURFACE 
FLAT: FLAT SEA SURFACE 

5-4 : SEA SURFACE ROUGHENED 3Y BEAUFORT 4 WIND 

Ts : CALCULATED AEPARENT SEA SURFACE TEMPERATURE 

Td : DIFFERENCE BETWEEN ACTUAL AIsD CALCULATED SEA TEMPERATURES 

Tm : AGA SEA TEMPERA.TURE MEASUPvE.MENT WITHOUT ATMOSPHERIC ATTENUATION 

Tmd : DIFFE?E:NCE BETWEEN ACTUAL AND AGA SEA TE.MPERATURES 



ALL TE.MPERATURES ARE IN ®C 



112 



WIND Sl'KKI) TO l)T RATIO KKLATION 
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I'igurc A. 2 AT ratios plotted versus wiiul speetl. 



APPENDIX B 

IN-BAND FLUX CALCULATION USING TI-59 PROGRAM 



1, GENERAL 

The program for the TI-59 calculator used to calculate the in-band flux makes 
use of the Simpson's approximation prograjm included in the master library module. It 
was originally written for the HP-67 calculator by Robert Pitlak ( Appolo Lasers, 6357 
Arizona Circle, Los Angeles CA., 90045 ). and has been translated for the TI-59. 

To use this program the following values have to be inserted: 

a. The initial and final wavelengths of the spectral band. 

b. The temperature of the object ( in °K or °F or °C ). • 

c. The even number of intervals to be used bv the Simpson's approximation of the 
integral. 

The calculated flux is given in Watts; cm^ or photons/cm~sec 



2. FORMULAE AND CONSTANTS 

As we discussed in Chapter VI, Section E, Planck's blackbody law can be written 
in the form 

Q = J Cq q dX photons'cnmsec 

>•1 

X. 

or XI = J Cj (q’X) dX Watts cm^ 



1 C, 

Where q = x = 

X^(e-^-I) 

Cq = 2kc = 1.883651556 X 10'^ ( pm'^/cm'sec ) 
Cj = 2rthc2 = 37418.42875 ( Watts ft ) 
C 2 = ch k = 14388.32334 ( pm °K ) 



The integrations are done numerically using Simpson's approximation. 
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3. PROCEDURE 



How to use the program 

a. Load the magnetic card 

b. Press RST 

c. Enter initial wavelength in }.im. press R S 

d. Enter /o final wavelength in jim. press R S 

e. Enter temperature (T). press 

(1) Aifin’K 

(2) Bifin^F 

(3) Cifin“C 

f Enter even number of integration intervals 
g. Press 

fl) D to get the answer in Watts cm" 

(2) E to get the answer in photons, cm~sec 

4. PROGRAM LISTING 

In the following Table 24 we give the listing of the program used in the Tl-59 
calculator. 
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TABLE 24 



TI-59 PROGR.‘\M LISTING 
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02 2 


1 1' 4 




lm:: 








055 


07 


7 


115 


03 3 


175 


•m ^ 


370 








054 


03 


0 


114 


03 3 


r 4 


04 


04 








057 


?5 


a 


11 : 


04 4 


1 


34 


r CM 








053 


42 


3TD 


113 


55 • 


173 


0? 


0? 


000 


42 


370 


05? 


00 


00 


11? 


43 FCL 


17? 


14 


D 


00 i 


01 


01 


040 


?1 


F/3 


120 


07 07 


130 




I MV 


002 


?? 


PPT 


04 I 


7 4 


LPL 
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5ci - 


1 3 1 
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13 
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005 
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02 
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00 


00 
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1\ 
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0 • j b 


■? *? 


?RT 


04 5 


01 


1 


125 


54 > 


135 


1 


L2L 


66 r 


91 


P.'3 


0 4 c* 


05 
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1 2 1' 


75 - 
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4 e 

i, >J 


£ 
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74 
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00 


0 
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01 1 
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?3 
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00’? 


1 i 


R 
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00 


0 
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54 ) 
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1 -Ti 
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0 1 ' j 
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00 


0 
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54 > 
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05 


05 
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00 


00 
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0 
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1?0 
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01 
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04 


04 
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00 
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023 


0 4 
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53 ( 


203 


42 


370 


625 


00 


00 


034 


03 


3 


144 


43 RCL 


' 204 


03 


03 


023 


?1 




035 


?5 


a 


145 


02 02 


205 


01 


1 


027 


7 4 
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